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jj^ , Abstract. The parametric equations of the surfaces on which highly resonant quasi- 

periodic motions develop (lower- dimensional tori) cannot be analytically continued, 

■ in general, in the perturbation parameter e, i.e. they are not analytic functions of 
e. However rather generally quasi-periodic motions whose frequencies satisfy only 
one rational relation ( "resonances of order 1 ") admit formal perturbation expansions 
in terms of a fractional power of e depending on the degeneration of the resonance. 

\^ ' We find conditions for this to happen, and in such a case we prove that the formal 

. expansion is convergent after suitable resummation. 
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o 1. Introduction 

| Resonances play an important role in the theory of dynamical systems. A possible application 

is provided by problems of celestial mechanics, such as the phenomenon of resonance locking 
between rotation and orbital periods of the satellites [GP]. In fact, the presence of friction can 
select resonant motions which remain stable when friction (on astronomical time scales) becomes 
negligible. In such a case maximal KAM tori can be really observed only approximately and on 
very short time scale, whereas, on very large time scales one expects that only periodic motions 
survive. On intermediate time scale one can imagine that quasi-periodic motions, involving a 
number of frequencies less than the total number of degrees of freedom (and decreasing with time) , 
describe most of the observed dynamics. This makes interesting and important to study quasi- 
periodic motions occurring on lower-dimensional tori for nearly-integrable Hamiltonian systems. 
These quasi-periodic motions are characterized by frequencies satisfying s rational relations, with 
r = N — s ranging between 1 (periodic motions) and the number N of degrees of freedom (KAM 
tori). The number s equals the number of normal frequencies appearing in the perturbed motions, 
whereas r is the number of independent components of the rotation vector. 

The analysis of such motions is simpler under some (generic) non-degeneracy assumptions on 
the perturbations. On the contrary the situation becomes immediately very complicated if no 
restriction at all is made on the perturbation. Situations of this kind arise also in similar contexts: 
we can mention the conservation of KAM tori under perturbations for systems of N harmonic 
oscillators, proved for N = 2 but conjectured to hold in general [H], [R2], and the stability of 
Hill's equation under quasi-periodic perturbations [GBC]. In the case of lower-dimensional tori, 
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the non-degeneracy assumption is that the normal frequencies become different from zero when 
the perturbation is switched on. If such an assumption is removed, then only partial results hold, 
and only in the case s = 1, that is only in the case of one normal frequency [Chi], [Ch2]. 

Starting from the work of Eliasson [E2], a new approach to KAM theory of quasi-integrable 
Hamiltonian systems arose, based on the analysis of cancellations in the "Lindstedt series" for the 
functions mapping unperturbed motions (uniform rotations, in suitable coordinates) into corre- 
sponding perturbed ones. 

A convenient way to exploit cancellations to resolve apparent divergences in the Lindstedt series is 
through methods inspired by quantum field theory, consisting in graphical expansions, summation 
of classes of diverging subdiagrams, iterative study of the flow of the effective constants and the 
possible introduction of counterterms. With these techniques and ideas, a number of known results 
have been reproduced and new ones have been obtained; see [Gl], [Ge2], [Ge3] and [GBG] for some 
reviews. 

In this paper we follow the latter approach to investigate the conservation of (N — l)-dimensional 
tori for systems with N degrees of freedom. More precisely we consider N degrees of freedom 
systems described by analytic Hamiltonians of the form 

H(I,<p)=H (I)+ef(I,<p), (1.1) 

with (I, ip) G V x T N , where T = R\27rZ is the standard torus, V is an open subset of M N , e is a real 
parameter, and the free Hamiltonian H (I) is assumed to be uniformly convex: d]H (I) > C > 
for all leV. 

Definition 1 (Simple resonance). A simple Diophantine resonance for the unperturbed system is 

a motion taking place on the torus {I } x T N with Lu d = diH (I ) satisfying a rational relation 
uiq ■ i/q = for some v$ G Z W and \uj$ ■ v\ > C\v\~ T for suitable C, t > and for all ^ v G Z W 
not parallel to vq . 

Thus, if e = 0, the motions with rotation velocity ujq will foliate the torus {/ } x T W into a one 
parameter family of invariant tori of dimension N — 1. 

For e^O invariant tori with dimension N — 1 run by quasi-periodic motions with spectrum w 
will, in general, only continue to exist "close" to some of the unperturbed tori {Iq} x T N . The 
problem is simpler under non-degeneracy assumptions on the average (/) of the perturbing function 
/ on the torus {/ } x T^; it has been studied in [GG1] and [GG2] with techniques employed, under 
different assumptions on the perturbation, in this paper. Define the average of / on {/ } x as 

(/(¥>)>=' Hm i f ' f(<p + u t)dt. (1.2) 

It depends nontrivially on in the sense that in general it is a non-constant periodic function of 
ip. This is more easily visualized in coordinates adapted to the resonance: they are defined by a 
linear canonical transformation (/, ip) d = [Iq + 5~ 1 (A, B), S T (a 7 /?)) , with S a non-singular integer 
components N x N matrix with determinant det S = 1 such that ojo = S T (u>, 0) with u> G IR^" 1 , 
\u-u\> C \v\- T « for allO v G and A G IR^ 1 ,^ G E,ae T N ~\(3 G T. 
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In the coordinates (A, B,a, (3) the Hamiltonian becomes an analytic function of (A,B), in the 
domain obtained from V under the transformation S, and (a, (3) G T w_1 x T. It will be of the form 
H(A, B, a, (3) = H' Q {A, B) + ef(a, 13), where H' (A, B) = u> • A + H (A, B), with H vanishing to 
second order at A = 0, B = and uniformly convex in the domain where it is defined. Of course 
the functions H,H ,f have a different meaning with respect to those appearing in (1.1), but we 
prefer to use the same notation for simplicity. For the same reason we still shall use the notation 
(/, ip) to denote the new action-angle variables, by setting / = (A, B) and ip — (a, (3). 

So in the new coordinates the Hamiltonian H, the rotation vector u> G R w_1 and the average 
/o(A, B, f3) of the perturbing energy can be supposed to be such that 



H = < J j-A + H {A,B)+ef{A,B,a,(3), \u ■ v\ > ^ VO^eZ^ 1 , 
fo{A,B,/3) = f f d " f(A,B,a,/3), if f(A,B,a,0)= £ e iu ' a f„(A, B, (3), 



(1.3) 



near the unperturbed resonance and without any loss of generality. The corresponding equations 
of motion for X(t) = (A(i), B(t), a(t), (3(t)) are 

A = — ed a f, 

B = -ed f, ■ /n -l\ 

A t, ^ n * ° r X = Ed x H(X), with E=r 1 , (1.4) 

a = u> + d A H + ed A f, V 1 u / 

$=d B H Q +ed B f 1 

where E is the standard 2N x 2N symplectic matrix. We study the existence of motions which 
can be described by a constant (3o and by functions A(ip),B(ip),a(ip),b(ip) of ip G f Ar_1 ; which 
tend to as e — > 0, such that, by posing 



A(i) = A(^ + wi), B(t) = B(ip + ut), 

a(f) = ^ + uf + a(t|) + uf), 0(t) =#> + b(^ + wi), 



(1.5) 



one obtains, for all if> G T^ -1 , solutions to the equations of motion (1.4). For brevity we shall 
sometimes write X(t), instead of X(ip + cot), to indicate solutions of (1.4) of the form (1.5). 
Note that, for e = 0, the motions (1.5) reduce to 

A< >(t) = 0, S(°)(t) = 0, a<®(t) = il> + ut, (3W(t)=(3 , (1.6) 

where ip G T^" 1 are arbitrary. The motions X^(t) = (A<°) (t), (t), a<°> (t), (3^ (*)) represent 
the unperturbed resonant motions filling a one parameter family of AT — 1 dimensional invariant 
tori (parameterized by /3q). 

Given any function G(ip, e) we shall denote by G v the i/-th Fourier component of its Fourier 
expansion in ip and we shall call [G] k the k-th. order term obtained by expanding in e the function 
G. Furthermore we shall denote by [G] k the i/-th Fourier component of the Fourier expansion of 
[G] k and by G^o the function G — Gq. 
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A formal solution of (1.4), (1.5), as a power series in e, X(ip) = + eX^{tp) + . . ., is 

well known to exist if f3 G T is a stationarity point for the average /o(0, 0,(3) (i.e. a point with 
such that 9^/0(0,0,^0) = 0) which is not degenerate (i.e. <9|/o(0, 0, /3 ) 7^ 0); cf. for instance [P], 
[JLZ] and [GG1]. Here we consider explicitly the case in which the non- degeneracy condition on 
<9|./o(0, 0, po) can fail t° hold. This is a case in which in general no formal solution in powers of 
e can be constructed. 

If (3q is such that dpfo(Po) = 0, then there exists a function X' such that X^ + eX' solves (1.4) 
up to terms of order e 2 (excluded). Namely X'{ip) = (A'(ip),B'(ip),a.'(ip),b'(il))) is obtained by 
applying the operator (u; • e?^,) -1 to the vector 

( -<V(O,0,V,A>) , 9 / /(0,0,V,/3o) + 9/ 2 ffo(0,0)(A'(V),S'(T/>))), (1.7) 

where d v ,di denote respectively the derivatives with respect to the angle and action variables. 
This means that one first determines A'(ip) and B'(ip) by solving (u> ■ d^,)(A'(ip), B'(ip)) = 
— 9^/(0, 0,t/>,/?o); the (otherwise arbitrary) averages A' Q ,B' Q are fixed by requiring that the op- 
erator (u> ■ 9) _1 can be applied to the vector formed by the last N components of (1.7), i.e. 
9 7 /o(0,0,/3 ) + 9|F (O,0)(A o ,B o ) = 0. In this way also a'(ip) and b'(ip) can be obtained. The 
averages of the angle variables will be chosen a = while we leave b' as a free parameter (to be 
suitably fixed at higher orders to make the equations (1.4) formally solvable). 
There are, however other solutions which are correct up to order e 2 (excluded). If 9j/o(0, 0, (3q) = 

for all j < ko and 9^ 0+1 /o(0, 0, /3q) 7^ 0, one can imagine to add to the free parameter b' any 

1 

polynomial in powers of r\ = |e| fc o of degree < ko- It can be checked (and it will be explicitly 
shown in next sections) that, for any choice of this polynomial, a solution of the equation of 
motions correct up to terms of order e 2 (excluded) exists. This is, ultimately, the reason why 
a consistent expansions in powers of e cannot in general be continued beyond first order: the 
consistency condition of the equations of motion necessary to improve (to second order included) 
the solution fixes one coefficient of such (a priori arbitrary) polynomials in r\ and, in general, this 
really forces the expansion to be an expansion in powers of r\ different from 77^° = e. 

The analysis below shows that an expansion in r\ is actually possible at all orders if the following 
assumptions are satisfied. 

Assumptions. 

(a) The constant a = J [dpf(XW + e X )] is a ^ 0. 

(b) The matrix <9fiJ (0,0) is positive definite. 

(c) There is k > such that 9^/ (0, 0, (3 ) = for all j < k a and c d = ^9^ 0+1 /o(0, 0, (3 ) + 0. 



The constant a appearing in assumption (a), written explicitly, is 

S Av /(I (0) (^)) • (« • ^)- 1 (9 7 /(x(°)(V)) +9 2 ff o (O,0)(A'(</0,B'W)) - 



def 

a = 



9fti/(^ (0) W) • (« • d^)- l d v f{X^{^)\ 



^° (1.6 





2 , 



For instance if / depends only on the angle variables tp and d 2 H(0, 0) = 1 the constant a is 

W\ 2 \UP)\ 2 + \d P UP)\ 2 



(1.9) 
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The number fc appearing in assumption (c) is a measure of the "degeneration" of the resonance, 
while the order of the resonance is the number of rational relations between the unperturbed 
frequencies, which is 1 in our case as the unperturbed motions have N — 1 independent frequencies. 
The case fco = 1 was considered in [GG2] and we will not consider it again here. We shall focus 
on the case fco > 2, in which case a formal power solution in e to (1.4) docs not exist. In fact, one 
checks that, as a consequence of assumptions (a) and (c), the average over xjj of the r.h.s. of the 
second equation in (1.4) is different from at order e 2 , for any possible choice of the free parameter 
b' introduced after (1.7). 

However, under the assumptions above, it is possible to find a formal solution to (1.4) such that 
the average of b is b = O(|e| 1//feo ): the average of b will be fixed in terms of the constants a and c 
in such a way that the average over tp of the r.h.s. of the second equation in (1.4) is at all orders 
in lei 1 /' 00 . 

The necessity of a fractional powers expansion can be seen by a heuristic argument sketched 
in Appendix Al: the argument also suggests, as a conjecture, the forthcoming Theorem 1 and 
motivates the assumptions (a) to (c). 

In Appendix Al it is in fact shown that, in the simple case Hq(A,B) = \{A 2 + B 2 ) and / 
depending only on the angles, a canonical transformation (explicitly constructed in Appendix Al), 
defined in a neighborhood of {0} x {0} x T^ -1 x {/3 } 7 maps w • A + H (A, B) + ef(a, (3) into 

w • A + H (A, B) + - f3 ) ko+1 +e 2 a(f3- f3 ) + 

fc o + l (1.10) 

+ 0(e(J3 - /?o) feo+2 ) + 0(e 2 (fi - (3 ) 2 ) + 0(el 2 ) + 0(s 3 ). 

As discussed in Appendix Al, the Hamiltonian equations corresponding to (1.10) can be consis- 
tently solved to order e 2 . In particular, for some choices of the signs of e, a, c, the angle (3 admits 
an approximate quadratic equilibrium point 0(|£r| 1 / fe ° ) , whose stability depends again on the rela- 
tive signs of e,a,c. This second order computation suggests the conjecture that the unperturbed 
motion X(°\t) can be continued at e ^ 0, provided the average of (3 is chosen O(\e\^ ko ). The 
perturbed motion, if existing beyond second order, will take place on a torus (a small perturbation 
of the free one) which we shall call elliptic or hyperbolic, depending on the stability of the behavior 
of the linearization of the motion of (3 in the vicinity of its equilibrium point: if the corresponding 
pair of nonzero Lyapunov exponents is imaginary then the torus will be said elliptic, if it is real it 
will be said hyperbolic. 

In the following a sparse Cantor set dense at will mean a set £ contained in an interval 
/ = [— £0)0] or I — [0, £o] with an open dense complement in / and with as a density point in 
the sense of Lebesgue integration. In particular £ will have positive measure. 

Theorem 1. Consider the system described by the Hamiltonian (1.3), under the assumptions (a), 
(b) and (c). There exists Sq > such that for \e\ < So the following holds. 

(i) If fco is odd and ce < there is at least one hyperbolic invariant torus of dimension N — 1 
with rotation vector lj. If ce > there is a sparse Cantor set £ C [— £o,0] with the property that 
for e G £ there is at least one elliptic torus of dimension (N — 1) with rotation vector lj. 

(ii) //fco is even and ec has the same sign of —a, there is at least one hyperbolic invariant torus 
of dimension N — 1 with rotation vector u). Moreover, there is a sparse Cantor set £ dense at 
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such that, if e G £ and ec has the same sign of —a, there is at least one elliptic invariant torus of 
dimension N — 1 with rotation vector uj. 

Remarks. (1) If assumption (c) is violated, since / is analytic in (3, then /o(0, 0, (3) = as a function 
of /?. In this case we can perform a canonical transformation removing the perturbation at order 
e and casting the Hamiltonian into the form H' = H' + e 2 f, for some new analytic functions 
Hq and /'. If /' satisfies the assumptions above we can apply Theorem 1 to H' . If H' satisfies 
assumption (c) and violates assumption (a) we cannot say much. If H' violates assumption (c) we 
can again remove the perturbation at lowest order through a new canonical transformation and 
cast the Hamiltonian into the new form H" = Hq + e 4 f" and hope to be able to apply Theorem 
1 to H". And so on. 

(2) Assumption (a) is essential and, if it does not hold, our expansion may fail to be convergent. In 
fact, in some cases (fco even), it is easy to show that if assumption (a) fails there cannot be perturbed 
motions of the form (1.5), see Appendix A2 for an example. If k is odd, [Chi] proved existence 
of hyperbolic tori even if assumption (a) fails; in these cases we expect that a new perturbation 
parameter must be identified. The heuristic analysis in Appendix Al concretely suggests plausible 
results to be expected if a = (under alternative assumptions). 

(3) We expect that assumption (b) is not essential and that it could be weakened into the request 
that d]Ho(0,0) is non-degenerate and that 9^ B i?o(0,0) ^ 0. Certainly the convexity assumption 
simplifies some of the estimates (see Appendix A6) and we did not attempt to eliminate it. 

(4) The only known result on the problems considered here is in [Chi], where conservation of 
(N — l)-dimensional "hyperbolic tori" is proved under weaker assumptions, although we study 
conservation of both hyperbolic and elliptic (N — l)-dimensional tori under assumptions (a) to (c) 
above. 

(5) In principle one could proceed in a different way rather than following our approach. One 
could first perform the canonical transformation described in Appendix Al and leading to the 
Hamiltonian (1.10), then study the system so obtained with other techniques, such as those in 
[El], [Po], [Rl] or [JLZ]. To this aim one should use that in the new coordinates the unperturbed 
Hamiltonian contains terms of order 77 2,So_1 , w hile the perturbation is of order r/ 2k ° , hence has a 
further rj. 

(6) The analysis via Lindstedt series and summations of classes of diagrams is the main aspect of 
this work. Also the analyticity properties in 77 at a; fixed are, to our knowledge, new. 

(7) Technically the present work is strongly inspired by [GG1], [Gel] and [GG2]. The proofs that 
can be taken literally from [GG2] will not be repeated here: hence familiarity with the latter 
reference is essential. 

(8) The resummed series has manifest holomorphy properties which show that the set £ is in the 
boundary of a complex domain where the functions X{ip) are analytic in e: we do not discuss the 
details (sec [GG1],[GG2]). 

The paper is organized as follows. In Sections 2 and 3 we prove formal solvability of the equation 
of motion in power series in lei 1 /* and we describe a graphical representation of the terms appearing 
in the formal power series. Such a representation involves labeled rooted trees, and is very similar 
to diagrammatic representations through Feynman diagrams arising in quantum field theory. In 
Sections 4 to 6 we describe an iterative rcsummation scheme which eliminates some classes of 
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divergent subdiagrams and iteratively changes the power expansion, and we prove convergence of 
the resummed series. Some details of the proofs are deferred to the Appendices. 

2. Lindstedt series 

We define e — an k °, with a e {±1}, i) > 0, and look for a family of formal solutions of the 
equations of motion in powers of r\ parameterized by ^ £ T _1 , in the special form X(t) = 
(A(t),B(t),a(t),(3(t))- 

OO 

A(t)= E ^ E e^'A? 1 . 

OO 

B{t) = E V h E e iv ^+^B^\ 

OO 

a(t)=V + «t + a(t), a(t) = E r? fe E e^'^+^a^ >, 

k—ko v^O 

OO OO 

0(t) = A, + 6(t), &(*) = E ^ E e^-^+^ftW + E »M" 

fc=/c i/^O fc=l 

where all the involved functions have also a dependence on (3q which has not been made explicit. 
The formal series (2.1) will be called Lindstedt series as it extends the corresponding notions 
already used in the theory of quasi-periodic motions on maximal tori. 

In the following the average will be abbreviated as 

Remarks. (1) The functions A,B,& and b — bo have been chosen as power series in n starting 
with the fco-th order as the first non trivial order; a has been chosen with zero average (this just 
corresponds to a redefinition of the origin of T W_1 ), while the average of b has been chosen as a 

series in r\ starting with the first order; the coefficients f3k d — b k ^ will be chosen in such a way that 
the Lindstedt series admits a formal solution. 

(2) With the choices in (2.1), equations (1.4) are identically solved for any order k < k in n. The 
parameters /3fc, fc < fc , are left as free parameters, to be explicitly chosen below. 

To write the generic fc-th order of (1.4) we introduce the following definitions: given any function 
tp — > F(X(ip)), let [F]( k ^ be the fc-th order in the Taylor expansion of F in n and let [F]^ be 
the f-th Fourier component of the Fourier expansion of [F] ^ . Note that this notation concerns 
expansions in n and is different from the one introduced after (1.6) which dealt with expansions 
in e: with the new notation the quantity denoted a = [dpf(X^) + e X')]q before (1.8) has to be 
written as a = a [dpf(XW + an k ° X')]^ as the superscript k denotes fc-th order in s while ( fe ) 
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denotes fc-th order in 77. Then, if v ^ and k > k , the equations (1.4) become 

(iu ■ u) a« = [8 A H Q ]^ + o-[d A ft- k °\ 
{iu ■ v) = [d B H l k) + <y[d B f] { t k °\ 
and, since uj satisfies the Diophantine property, they can be solved provided 

0=[dp.f]t ka \ (23) 

0=[d A H i k) +cr[d A fi k - ko \ 

0=[d B H i k) +a[d B .fi k - k « ) . 

If assumptions (a), (b) and (c) in Section 1 are satisfied, such a formal solution can be shown to 
exist provided the formal series for the average 60 is suitably chosen. 

Lemma 1. Under the assumptions (a), (b) and (c) if ko is odd, a formal solution of (1-4) in the 
form (2.1) always exists; if ko is even, a formal solution of (1-4) in the form (2.1) exists ifaac < 0. 
When a formal solution exists, A^ k \ B^ k \ a^ k \ b^ are uniquely fixed in the case ko odd. If ko is 
even, there are two possible such sequences corresponding to the choices fii = b^ = ±(— aa/c) 1 ^ . 

Proof. Let X(°\ip) = (0, 0, ip, fio), e — ar/ k ° with a = sign(e) and n > 0, and look for a formal 
solution t — > (A(t), B(t),at(t), (5{t)) obtained by setting ip — rp + ujt in 

00 00 

k—ko k—1 

where ^ = (0, 0, 0, (3 k ). Set = for < h< k and = 0. 

Suppose inductively that X^> has been determined for h = 0, 1, . . . , k — 1, for k > 1. Consider a 
generic polynomial Y(ip) d = Ylh=o n h Y^ h \tp) and suppose, inductively, that if Y^> = X^ + 
the function Y(u)t) solves the equation of motion X = EdxH(X) up to order k — 1. This is true 

for k = 1. Let A d = (w • 9,/,) and remark that the identities fY^j=i ^tpi^j ' &YjH(Y) dip = and 
/ Sj=i ^V>« ^j' ' dtp = arc identities /or a/Z periodic functions Y(tp). 

Then = d i>i Y i ' (( S AY )i + d Y 3 H{Y)) dip and Y = Ed Y H(Y) + 0(r] k ) imply (note that 
Y has degree < k in »;) / d^*^ • [<9y 3 #(y)]( fe ) dip = 0, i.e. /[0 a JT(y)]( fc ) dip = or 



fe— 1 

[d a f(J2 V h (XW +^)))]( fe - fe °) =0 , (2.5) 

/t=0 

which is the first of the compatibility conditions (2.3). This leaves a^ as an arbitrary parameter: 
we can set a k ^ = 0. We also remark that, if I = (A, B), the last two of the four equations in (2.3) 
have the form 

{Ed) v d I H {0,Q)I^ + F^({x( h \^} <h<k-i) = 0. (2.6) 
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By assumption (b) in Section 1, d]H (O,0) = (Ed) v diHo(O,0) is invertible; hence to impose 
the last two compatibility conditions in (2.3), it suffices to fix I^ k) = -[d]H o (O,0)}- 1 F < - k '>({X ( - h \ 
£^}o<h<fc-i)j which is a known function of {X^ h \ ^ h ^} a <h<k-i, for any possible choice of the 
functions {£ ( ' l) }o<h<fc-i- 
So, to fulfill the compatibility conditions (2.3) and to continue the inductive construction of 
we are left with imposing the second of (2.3), which can only hold if Pk-2k a +i satisfies 
certain compatibility properties. In fact if k < 2fco there is no requirement because the necessary 
condition that [dpf] < £~ ko * > = is automatically satisfied by assumption (c). 

For k = 2fc the condition [dpf} [ Q 0) = can be expressed as follows. Let c = -^dp 0+1 fo(0, 0, 
and remark that with the notations leading to (1.8) it is X^ ka ^ = aX': then the second of (2.3) 
becomes 

d px f(X^)-X^ + ^^ +1 /(I (0) )(ft) fc »| = aa + cf3 k ° = 0, (2.7) 



fco! 



which means 



q _ f {-ao-/c) 1 / ko if k is odd, ^ ^ 

' \ ±{-aa/c) l / ka if k is even and aac < 0, { ' ' 

which is the compatibility condition to which (3i must be subject. 

For k > 2fc the second of (2.3) only involves a sum of quantities depending on X^ with 
h < k — k a and on with h < k — 2fc with the exception of a single term, proportional to 
-2fc +i) involving ^ k 2fco+1 ). Therefore the second compatibility condition in (2.3) can 
be fulfilled by properly fixing /3fc-2fe +i i n terms of X^ with h < k — fco and of with h < k — 2ko, 
provided (3\ exists and (3\ ^ 0, i.e. provided a ^ as assumed here. 

This means that if fco is odd all (3k arc uniquely determined while if fco is even there are two 
possible sequences /3k depending on the two choices for (3\ in (2.8). ■ 

Remarks. (1) For fc even the choice of j3k will be possible only if e has sign a such that —aa/c > 0. 
(2) With the notation (2.1) one has X {k) = (A fe) , B {k) , 0, 0) as a fc) = and b {k) = ^ k \ 



3. Tree formalism and formal series 

Given that the formal Lindstedt scries is well defined, by proceeding as in [G2] a graphical repre- 
sentation of the contributions to the perturbative series will be introduced. The idea to get the 
rules explained below is to start by representing the r.h.s. of (2.2) as a power series in r\ by simply 
expanding the functions d v f,djH in their arguments around X^ (ip) and then each argument 
again in powers until one obtains a power series in rj. 

Since the arguments of / depend on the X^ with h < k we obtain recursively a natural 
representation in terms of trees. Consider (X^) 7 where 7 ranges in the symbols list 

def 

A = {Ax,.. . ,A N -x,B,ai, . . .,a N -!,(3} . (3.1) 



This notation is more convenient than the alternative one which would simply label the components 
with a label 7 = 1,2, ... ,N, N+l, . . . , 2N because the first N — 1 components play a very different 
role than the N-th or the remaining ones. It can be disturbing as the labels in (3.1) have the 
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meaning of canonical variables when they do not appear as labels: hence the reader should keep 
this in mind in what follows. Occasionally, only in cases of possible confusion, we shall use the 
labels 1, . . . , 2N instead of (3.1). 

It will be convenient to call the first N labels action components and the last ./V angle components; 
the first N — 1 of the two groups will be called, respectively, fast actions and fast angles component 
labels while the last will be called a slow action or, respectively, slow angle component label. When, 
occasionally, it will turn out to be useful, we shall denote the N action labels simply by I and the 
angle labels by <p. If 7 G A and 7 is one of the first N labels (i.e. it is an action component), then 
7 + N will indicate the corresponding label of the angle components and vice versa if 7 G A is an 
angle component then j — N will denote the corresponding label in the group of action components. 
With a slight abuse of notation, given 7 G A, we shall write 7 G / if 7 is an action component and 
7 G p if it is an angle component. If 7, 7' G A then <5 7 y denotes the Kronecker delta; we shall also 
use the notation 

6 ™ = {o, if7G/,' 6 ^ I = {o, if^e/. (3 ' 2) 

Looking at (2.2) we see that, for v ^ and k > ho, we can write it as 




Here the symbols (9 7 have to be interpreted as derivatives of f v with respect to action arguments of 
/ if 7 is an action component, as derivatives with respect to f3 if 7 = [3 and as multiplications by ivi 
if 7 = a>i, i = 1, ... ,N — 1. A summation convention is conveniently adopted for pairs of repeated 
component labels. In the second line no £( fe ™) appears because h as 0IUV an gl e components but 
H does not depend on the angles. 

For v = and 7 G / (i.e. 7 an action component; cf. Remark (2) at the end of Section 2), we 
use (2.6) and write 



n=l 



00 p 

(4 fc) ) 7 = - (d*H (0, 0)- 1 )^,.^ [£ £ ^(250)y0 7l ... 7j> ff o (O,O) \{{X^) ln + 

(3.4) 



p—2 ki + ...k p = k 
!>! + .. . + l> v =0 



ff E E 3(^) r 9 7l ... 7p /, (0,0,/3o)n(^ n) +^ fen) ^no) 7n 

p— 1 fci + .-.fc p = fc-fco n— 1 

^0 + ^l + --- + ^p=° 



For z/ = and 7 an angle label the (2.2) are identically satisfied if the are determined as 
prescribed in Section 2. The equation that fixes is again a relation of the type of (3.3) and 
(3.4) as we shall discuss in detail later. 
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Proceeding as in [G2] and [GG1] we represent {xi h) ) 1 as 7 4^ • W and ^ = (0,0,0, 

as 7 ^ O ^ and we realize that (3.3) and (3.4) can be very conveniently (as it turns 

out) represented by graphs of the type represented in Figure 1. 




(kp+q) 



Figure 1. Graphical representation of equations (3.3) and (3.4). For i/^o the graph 
represents (3.3), while for i/=o it represents (3.4). 

The label <5o = 0,1, that we call the degree label of the node to which it is associated, identifies 
the terms in (3.3) or (3.4) which contain inside the square brackets H (then So = 0) or / (then 
Sq = 1). The other label attached to the first node indicates the harmonic Vo (mode label) selected 
in the terms with / (i.e. with 8q = 1) and we take Vq = if So = (because H does not depend 
on the angles). A further component label 7' is attached to the right extreme of the line exiting 
the central node in Figure 1 and it indicates the derivative {Ed)y in (3.3) and (3.4). We call a 
component label a "right" or "left" component label if it is attached at the beginning or at the 
end of the (oriented) line. 

The component labels attached to the first node determine the components of the tensors defined 

by the derivatives of /„„ or of Ho- The labels attached to the left extreme of each endline (to the 

right of the bifurcation point in the figure) determine which component of {X^"^) ln is taken in 

the products of X's in (3.3) or (3.4). Finally the root line symbolizes the factor that is outside 

the square brackets in the expression (3.3) or in the expression (3.4): it will be called propagator 
v 

of the line 4- — . Hence the propagator of the line will be the matrix 

(iuj ■ <5 7 y, if v 7^ 0, 



frr'M -{ -(d* Ho {o,o)\ . w m^> = (3 ' 5) 



and u will be called the momentum of the line. 

Finally the endpoints are divided into endpoints representing for some h, marked by bullets, 
and others representing marked by white disks that we shall call leaves. 

It follows from the analysis of Section 2 that, setting (3\ equal to the constant in (2.8) and 
replacing ( fco j_ 1 ^ ! ^ 0+1 /o(0, 0-/3 u ) with ko c by the definition of c before (2.7), the coefficients fih 
with h > 2 can be derived in terms of f3\ from the relation 

&-2fc„ + i = ~ L-x E E* ^(Ed) B d 11 ... 1 J Uo n (*& m) +i {km) 5u m0 ) lm , (3.6) 

&0 c Pi s >! fc 1 + ... + fc s = fc-fc ' m=l 

v Q + v 1 + ...+is s =0 
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where the derivatives of f Vo are, as above, evaluated at (0, 0, /3o) and the * on the sum recalls that 
we are excluding from the sum the contribution equal to — fc c /3^ 0_1 /3k-2k +i- 

At this point we can repeat the construction and replace each endpoint of Figure 1 ending in a 
node or in a leaf with a node into which merge several lines each of which comes out of a node 
or a leaf with some labels (h),u, 7 representing (xi/^) 7 for some h, v, 7 or, in the case of leaves, 

Using also the relations (3.4) and (3.6), the construction can be iterated until we are left with 
a tree 9 whose endpoints cither carry a degree label 1 or are leaves representing j3\: we denote 
by L(9) the set of endpoints representing such leaves. The constraint (3.6) can be automatically 
implied by imagining that the propagator of a line i with momentum u is 

' (iuj ■ u) -1 <5 77 <, if u ^ 0, 

^<H = r 77 " ' s 5 *» = »,**m, (3 . 7) 

-(djH (0, 0)" 1 ) , N 5 7l S Yv + "f/Vi ' * " = 0, * £ L{B), 



where is the node preceding I on r. We can write {X v ) 1 as a formal power series in 77, whose 
terms can be computed in terms of tree values. Let a tree 9 be a tree diagram with nodes which 
look like the node drawn in Figure 1 and let V(9), L(9) and A(0) denote the sets of nodes, leaves 
and lines of 9, respectively. The tree value Val(0) will be a monomial in 77 obtained by multiplying 

(1) a factor 

[F(»«)]y e , 7 (0) =vr 1 ko (Ed)y e 9 7ol ... 70po /„ (0,0,/3o), if S = 1, 

[^)]^ ,7(») =(^ o 9 7oi ... 7OPB ^ (O,0), if<5 n =0, 

per each node o € mto which merge p > 1 lines carrying component labels 7(0) = 

(701 ■ ■ • 7op ) an d emerges a line i<o carrying a label 7^ ; 

(2) a factor an k °{Ed) 1 ' f Vv (0, 0, (3q) per each endnode D L(0) (note that necessarily <5 = 1); 

(3) a factor 77/^1 per each leaf o e L(6) (note that necessarily i/ = and jg B = B); 

(4) a factor g ln '(v) per line £ € A(0), given by (3.7). 

Note that the construction described above and in Section 2 forbids the presence in the tree 
diagrams of some configurations of nodes. More precisely, calling trivial the nodes t) with p v = 1 
and v v = and b-trivial the nodes c with p v = k , f = 0, 7^ = B and immediately preceded 
by at least fco — 1 leaves, the following configurations of nodes: 

(i) trivial nodes with (> = and the entering line with momentum and 

(ii) b-trivial nodes with the exiting line with momentum and all the entering lines with left 
component labels equal to (3 

are not allowed, in the sense that Val(0) = if 6 contains such configurations of nodes. The reason 
why the trees containing such configurations of nodes are forbidden is a consequence of the use of 
the relations (3.4) and (3.6). Trees with no forbidden node will be called allowed trees. 

The result is that the perturbed motion runs on a (n - l)-dimensional torus whose equations 
are formally written as a sum of values of (allowed) tree diagrams, computable by using rules very 
similar to those listed in [GG1], with the value of an allowed tree 9 defined as 

v&m = it n e*)w m i n *)( n g *)> ^ 
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def 

where Gii = g ltl i(yt) with v t = = X)ro<o v k>i while F is the node factor (a tensor) defined 
in items (1) and (2) above. In the product all the labels are summed over, except for the root 
label 7^ , where £q is the line entering the root. We call 0£ v the set of trees with degree 
+ fc o Eoev(0) ^ - ( 2fc o - 1) J2ieA(0) ^Mje0 s Y e B = k, v ta = v and j ta = 7 (these are trees 
whose value is proportional to w k ). As in [GG2], we denote by Ok,v,f the set of trees with k nodes, 
and with labels u^ = v and je = 7 associated with the root line. 

The definitions above are given so that the formal scries for A 7 is given by the sum 



The above rules give, order by order in powers of n, the solution of the perturbed equations of 
motion. In particular one can check that the trees contributing a monomial in n of degree k > 1 to 
the conjugating function have a number of lines that is bounded above and below proportionally 
to k. This is a property extensively used in the convergence analysis, for instance to show that the 
number of non-numbered trees of degree k is bounded by a constant to the power k; cf. [GG2] for 
details. An explicit bound is > k/k n and < 3fcofc, see Appendix A3. 

4. Elimination of the trivial nodes 

The problem of proving convergence of the series just defined is very similar to that treated in 
[GG2] . As in [GG2] the difficulty is that even exploiting the cancellations analogous to those of the 
maximal tori case, we are left with tree graphs containing chains of subdiagrams with one entering 
and one exiting lines, carrying the same momentum, that we call "self-energy (sub)diagrams". 
Naively such subdiagrams are the source of bad bounds on the fc-th the contribution proportional 
to i] k to the series. Proceeding as in in [GG2] we iteratively resum such chains into "renormalized 
propagators" and change step by step the structure of the perturbation series. At each step we 
define different rules to compute the tree values: we assign at each line I a scale label [tig], with 
ri£ > —1, depending on the size of its propagators, and, at the n-th step (n = 0, 1, . . .), we will not 
allow trees containing chains of self-energy diagrams on scale > n — 1; at the same time we will 
assign to each line a propagator different from that in (3.7), depending on the value of its scale 
label, as explained below. 
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Figure 2. Examples of lowest order self-energy clusters on scale [—1]. In (a) one 
has 7 = B and 7' € 7 or vice versa, while in (b) one has 7 = B and 7' = (3. The 
order label <5 in the first graph can be 5 = 0, 1. The leaves in the second graph 
represent factors (i\ so that they contribute to the order of the self-energy clustser 
a power r) k ° ~ 1 . 



The first step is the removal of the trivial nodes and of a class of few other subgraphs that can 
be present in allowed trees, see Section 3, whose value is not 0. Assign, in a way similar to that 
described in [GG2], a scale label [—1] to the 0-momentum propagator and let S^_ 1 be the set of 
self-energy clusters T on scale [—1] and degree k (i.e. such that k = \L(T) \ + J2vev(T) ^» ^0 — 
(2fco — 1) X^eA(T) &viO&nP&"t' l B)i consistently with the notion of degree defined after (3.9). 

A cluster on scale [—1] is either a single node or a maximal connected set of nodes and lines, such 
that all the lines are on scale [—1]; the lines of non-zero momentum that enter or exit the nodes 
of the cluster are called external lines. A self-energy cluster on scale [—1] is defined as a cluster 
on scale [—1] such that there are only two external lines (one exiting and one entering), and they 
are connected to the same node of the cluster. This means that a self-energy cluster is formed by 
a node with 0-mode together with all lines (if any) of momentum 0, hence on scale [—1], which 
are linked to by a path consisting of lines of momentum 0; sec Figure 2. A cluster on scale [—1] 
with only one exiting and one entering line, which are not connected to the same node, will not be 
considered a self-energy cluster on scale [—1]. 

The self-energy value, i.e. the value Vt(ii) of the self-energy cluster T on scale [—1], is defined, 
similarly to (3.9), as 

^)=]a(M n ^)(^) il(t)i ( n *•*){ n g *)> ^ 

1 v oev(T) vev(T) oeA(T) 

where V(T), L(T) and A(T) are the set of nodes, leaves and lines, respectively, contained inside 
T. Then we can define 

00 

M l0] (v)=Y, E V t(v), (4-2) 

k=k T£5f_ 1 

where the sum runs over all the self-energy clusters on scale [—1]. 

Remarks. (1) The clusters T in (4.1) and (4.2) with k = fc correspond to the trivial nodes 
appearing in the tree expansion of previous section. In the next section the general notion of 
cluster will be introduced. 
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(2) The clusters on scale [—1] must contain only lines on scale [—1] and there are infinitely many 
of them and (4.2) can be illustrated by the two clusters in Figure 2. The self-energy values Vt{v) 
corresponding to the graphs in Figure 2 give the lowest order(in powers of rf) contributions to 
different entries of the matrix .A/!' ', as discussed in the caption of Figure 2. 

(3) Unlike the case in [GG2], the equation (4.1), defining the matrix M^(rj), is really an infinite 
series; but it is still convergent, if r] is sufficiently small. Convergence of the series defining M^ 1 
is a straightforward consequence of the fact that a self-energy cluster of degree k has a number of 
lines bounded by constant times k, see Appendix A3, and of the fact that the propagators of the 
lines with 0- momentum can be bounded by an O(l) constant if the line is preceded by a leaf or if 
7 = ij,7' = tfj and can be bounded by an 0(1) constant times rj~ 2ko+1 if the line is not preceded 
by a leaf and 7 = [3, 7' = B. 



^ [0] W=E E V T =( Q p _%) (4.3) 



By construction, M\^{rj) is real and has the following special structure: 

Q R 

P -Q t 
k=k res^_ 1 x 7 

where the superscript t denotes Hermitian conjugation and P, Q, R are N x N matrices which to 
lowest order in r\ have the form (with a natural meaning of the symbols, in agreement with our 
convention on the component labels (3.1)) 

P = tiH 0+ e$f Ot Q=(_ ea ° j/o ). *=(S -f^fVo/^- 1 )' (44) 

where fo has (0, 0, (3q) as arguments. The complete expression, including the higher orders in r], 
simply replaces the non-zero terms in (4.4) by convergent series that we can write 

P = 8lH 0+ eM II , Q=(_jy, * =(2 _^o°imJ' ( 4 - 5 ) 

where the entries M 7 y (77) are bounded uniformly in 77, for 77 small enough. They can be computed 
order by order using the rules of the last section. The vanishing entries of the matrix M^(n) 
remain zero to all orders: a property which simply follows from the definition of value of a self- 
energy cluster and from the remark that a derivative <9 7 with 7 = 0; acting on /o must be interpreted 
as a multiplication by 0. Note that, since P, Q, R are real and P, R are symmetric, (77) satisfies 
the following symmetry properties: 

EM l0] (r])E= [M m (f])] T , [M l0] (f])Y =M m {r ] ), (4-6) 

where * denotes complex conjugation and T transposition. A consequence of (4.6) is that M^E 
is Hermitian. 

Now, we formally resum the chains of self-energy clusters on scale [—1] into the new propagator 

^"J-s^iW (4 ' 7) 
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where x ^ 0. This means that we modify the tree expansion described in previous section: the new 
expansion will involve only trees not containing self- energy graphs on scale [—1] (and in particular 
containing neither trivial nor 6-trivial nodes) and with the propagator of the lines with non-zero 
momentum replaced by (4.7). It is not clear at all that the new resummed series is well defined: on 
the contrary it will become clear that it is affected by divergence problems similar to those of the 
original series. In some sense, we just eliminated a few of the possible source of problems (that are 
actually an infinite class of divergent sub-diagrams). However, the idea is to begin by eliminating 
this first few sources of problems and then, step by step, iteratively eliminate one after the other 
all possibile sources of problems (first the less "dangerous" and then the more and more dangerous 
ones). 

Certainly we must at least suppose that ix — M. ^ (rj) can be inverted: otherwise the values of the 
trees representing the new series might even be meaningless! To give a meaning to (ix — M^^tj))^ 1 
it is sufficient to impose det(ix — M\°\rj)) ^ for x ^ 0, by eliminating a denumerable dense set 
of values of rj. One can compute the determinant of (ix — A4^ (rj)), finding dct(ix — M.^(rf))) = 
-(ix) 2{N -^ ■ x 2 +erj k °- 1 M /3 { 3 (d BB H + eM BB ) + e 2 \M pB \ 2 , so that the condition of invertibility 
of (ix — M^(rj)) for x ^ becomes 

x 2 +er t ka ~ 1 M f} p(d 2 BB H + eM BB ) + e 2 \M f}B \ 2 ^0 (4.8) 

for all x = uj ■ v, v £ Z^ -1 . In computing det(ia; — M^(rj)) the property M p B — —M B p e M. has 
been used. 

If rj is chosen according to (4.8), we have that the norm of g'- ' (x; rj) is equal to the eigenvalue of 
[ix - M [°1 (rj)] E with smallest absolute value: this is because \\ix - M^(rj)\\ = \\(ix - M^(rj))E\\ 
and, as remarked after (4.6), (ix — M^(rj))E is Hcrmitian. Here and in the following we use the 
uniform norm: given a 2N x 2N complex matrix g, we define ||g|| = sup z6C 2N ) | z | =1 \g z\, where 

M = E i= i N- 

An approximate computation of the eigenvalues of (ix — M^(rj))E, see Appendix A4, implies 
the following result. 

Lemma 2. There exists an 0(1) constant p > such that if x 2 > prj 2k °~ 1 then the resummed 
propagator g^-°\x\rj) in (4-7) can be bounded as: 

g^\x;rj) ^maxj^f,^}, (4.9) 
k x Mi 

where ^ and p$ are, respectively, the minimum and maximum eigenvalues of P, see (4-3) and 
(4-5). 

From now on we shall proceed following [GG2]. First of all we assume that \e\ is in an interval 
(e/4,e] such that, by setting e = rf° , we can define the integer no through 

C ,2 2 -2(„ +l) < pI1 jk -l < C 2 2 -2no 5 ( 4 jq) 

with p as in Lemma 2. 
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In the first range of scales (in which x 2 > 2Cq2~ 2 " ) the small denominators can be bounded 
by the "classical" small divisor x 2 and we proceed as described in the section "Non-resonant 
resummations" of [GG2]. 

For smaller scales we shall see below that the small divisor will bounded below by an O(l) 
constant times min{x 2 , \x 2 + er] k "^ 1 \^ n \x; e)\}, with \^(x;e) a suitable O(l) function. How- 
ever here the distinction between the hyperbolic case (er] ko ~ 1 \^ (x; e) > 0) and the elliptic one 
(er] k °~ 1 \^(x;s) < 0) has to be made: in the hyperbolic case the small divisors will be always 
bounded by an 0(1) constant times x 2 , even for x 2 < O(ei] ko ^ 1 ), while in the elliptic case we 
will have to proceed differently, essentially as described in the section "Infrared resummations" of 
[GG2]. 

5. Multiscale analysis and non resonant resummation 

The resummations will be defined via trees with no self-energy clusters on scale [—1] and with lines 
bearing further labels. Moreover the definition of propagator will be changed, hence the values of 
the trees will be different from the ones in Section 3: they are constructed recursively as follows. 

We introduce a multiscale decomposition (see [GG2]): we call tp(D) a C°° non-decreasing compact 
support function defined for D > 0, 

V>(£>) = 1, for D > C 2 , ip(D) = 0, for D < Cg/4, (5.1) 

where Co is the Diophantine constant of u>, and let x{D) = l — tp(D). Define also ip n (D) — ip(2 2n D) 
and Xn{D) = x(2 2n D) for all n > 0. Hence ipo = 4>, Xo = X an( i 

1 = yj n (D(x)) + Xn(D(x)), for all n > 0, (5.2) 

for all choices of the function D(x) > 0: in particular for D(x) = x 2 , that we shall now use. 

A simple way to represent the value of a tree as sum of many terms is to make use of the identity 
in (5.2). The resummed propagator g^-^(x;w) d = (ix — Ad^^rj))^ 1 of each line with non-zero 
momentum (hence with x 0) is written as 

g^°\x; V ) = M* 2 ) 9 l ~° ] (x; V ) + Xo(x 2 ) g^°\x; V ) ^ g M {x - v ) + g {>^} ( x; v ), (5.3) 

and we note that g^(x; rj) vanishes if x 2 is smaller than (Co/2) 2 . 

If we replace each t/[-°] (x; n) with the sum in (5.3) then the value of each tree with k nodes is 
split as a sum of up to 2 fe terms which can be identified by affixing on each line with momentum 
i//0a label [0] or { > l}. Further splittings of the tree values can be achieved as follows. 

Definition 2 (Propagators). Let n n be an integer, and for 1 < p < n n , give 2N x 2N matrices 
M^(x;rj) satisfying the symmetry properties 

EM [p] (x;n)E = [M [p] (-x;n)] T , [M [p \x; rf)] * = M [p] (-x; n). (5.4) 
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Let M [n] {x;n) = M [0] (n) and M [ ^ n] (x; rj) = £™ =0 M [p] (x; rf). Define for n > n > 1 the propa- 
gators 

gin-Ufa n) d - f ^™ (a;2) Xm{x 2 ) 



ix - M^ n -^(x;r]) 

TTl-l 



{>„} ( } */ _[^x^J_ (55) 

— All-" ^(xjn) 

5 1 ' W ia:-MK»](a;;»7)' 
and <?'°'(x; n) = ^(z 2 ) {ix — M^iri))^ 1 . We call the labels [n], {> n}, [> n] scale labels. 

Remarks. (1) The matrices M^{x;n) will be defined recursively under the requirement that the 
functions (A, B, a, b) defining the parametric equations (2.1) of the invariant torus will be expressed 
in terms of trees whose lines carry scale labels indicating that their values are computed with the 
propagators in (5.5). 

(2) To have the propagators in (5.5) well defined, we have to eliminate for each value of p a 
denumerable set of e's, by imposing that {ix — M\- P \x; rj)) be invertible, in analogy with (4.8). 

(3) So far no can be any integer number. It will be fixed as prescribed after Lemma 2, in such a 
way that a bound like (4.9) will hold for all propagators with p < no- 

To define recursively the matrices we introduce the notions of clusters and of self-energy clusters 
of a tree whose lines and nodes carry the same labels introduced so far and in addition each line 
carries a scale label which can be either [— 1], if the momentum of the line is zero, or [p], with 
< p < no, or [> no], with n the same integer appearing in the statement of Definition 2 (still 
to be suitably fixed). Given a tree 9 decorated in this way we give the following definition, for 
n < no- 

Definition 3 (Clusters), (i) A cluster T on scale [n], with < n, is a maximal set of nodes and 

lines connecting them with propagators on scales \p], p < n, one of which, at least, on scale exactly 

[n]. We denote with V{T), L{T) and A(T) the set of nodes, the set of leaves and the set of lines, 

respectively, contained in T. The number of nodes in T will be denoted by kr- 

(ii) The itit > lines entering the cluster T and the possible line coming out of it (unique if 

existing at all) are called the external lines of the cluster T. 

(Hi) Given a cluster T on scale [n], we shall call nx = n its scale. 

Remarks. (1) The clusters on scale [—1] were defined before (4.1): they can contain either only 
lines with scale [—1] or no line at all (i.e. they can contain just a single node). 

(2) Here n < n . However the definition above is given in such a way that it will extend unchanged 
when also scales equal to or larger than no will be introduced. 

(3) The clusters of a tree can be regarded as sets of lines hierarchically ordered by inclusion and 
have hierarchically ordered scales. 

(4) A cluster T is not a tree (in our sense); however we can uniquely associate a tree with it by 
adding the entering and the exiting lines and by imagining that the lower extreme of the exiting 
line is the root and that the highest extremes of the entering lines are nodes carrying a mode label 
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equal to the momentum flowing into them (cf. [GG2], Figure 3). 

Definition 4 (Self-energy clusters). A self-energy cluster on scale [n], with n> 0, of a tree 9 will 
be any cluster T on scale [n] with the following properties: 

(1) T has only one entering line and one exiting line l\; 
( u ) Eoev(T) u * = 0; 

(iii) there is no line on scale [—1] along the path connecting £^ to i\. 
We call kx the number of nodes in V(T). 

Remarks. (1) The essential property of a self-energy cluster is that it has necessarily just one 
entering line and one exiting line, and both have equal momentum (because X^oev(T) u o = 0)- 
Note that both scales of the external lines of a self-energy cluster T are strictly larger than the 
scale of T regarded as a cluster, but they can be different from each other by just one unit. 

(2) The self-energy clusters on scale [—1] were defined before (4.1). 

(3) For n > 0, the number of nodes of any self-energy cluster on scale [n] is > 2, and the corre- 
sponding degree is > 2k . This can be seen as follows. Call T the connected subset of T containing 
no line on scale [—1] and containing the two nodes to which the external lines are attached. Then 
To must have at least two nodes with 5=1, and T \T is the union of subtrees with positive 
degree. 

(4) The clusters which satisfy properties (i) and (ii), but not (iii) are not considered self-energy 
clusters. The same happened for the self-energy clusters on scale [—1]. The reason for such a def- 
inition is that the cancellation mechanisms that will imply the bounds needed on the derivatives 
of the self-energy values (see next definition and Lemma 3 below) can be derived only under such 
an extra condition. On the other hand the clusters which verify only the first two properties, but 
contain lines on scale [—1] along the path connecting the external lines, require no resummation, 
and can be dealt with in the same way as the other clusters which are not self-energy clusters. We 
note that in [GG2] property (iii) was explicitly required only for self-energy clusters on scale [—1] 
but it was not mentioned any more for the others: the proofs in [GG2], however, implicitly used 
property (iii) also for the self-energy clusters on scale > — 1. 

Definition 5 (Renormalized trees). Let 0j^„ 7 be the set of trees with degree k (see comments 
after (3.9)), root line momentum u and root label 7 which contain no self-energy clusters. Such 
trees will be called renormalized trees. 

Definition 6 (Self-energy matrices), (i) We denote with S^ n the set of self-energy clusters with 
degree k and scale [n] which do not contain other self-energy clusters; we call them renormalized 
self-energy clusters on scale n. 

(ii) Given a self-energy cluster T e S^ n we shall define the self-energy value of T as the matrix 1 

^•^)=ia(W( n ^)(^) il(t)i ( n *)( n ^ ] ) ^ 



This is a matrix because the self-energy cluster inherits the labels 7, 7' attached to the endnode of the entering 
line and to the initial node of the exiting line. 
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where </|"*' — g^(u) ■ vf,rf). Note that, necessarily, ng < n. The fcy — 1 lines of the self-energy 
cluster T will be imagined as distinct and to carry a number label ranging in {1, . . . , kx — 1}- 
(Hi) The self-energy matrices M^(x;n), n>\, will be defined recursively as 

n—l oo n— 1 

MW(x;r)) = (l[x P (x 2 )) £ £ V T (*; 77) =' ( [] X P (x 2 )) 77), (5.7) 

p=0 k=2k a TeS^ n _ 1 P=0 

where the self-energy values are evaluated by means of the propagators on scales [p], with p = 
-1,0,..., n. 

The definition (5.7) makes sense because we have already defined the propagators on scale [0] 
and the matrices (x; 77) = M\ ^ (77) (cf. Definition 2). Of course we have still to check that the 
series converges. 

With the above new definitions let h vn with 7 = A, B, a, (5 be the values of A„, B v , a v ,b u . We 
have the formal identities 

00 

Vt = E E Val ^)' ( 5 - 8 ) 

where we have redefined the value of a tree 6 0^ as 

Yam= \Am( n ^)(^i) il(t)i ( n 9™io>"v,r,j){ n (5 - 9) 

1 ^ x>ev(B) £eA(0) oev(0) 

with [ne] = [— 1], [0], . . . , [no — 1], [> no]. Note that (5.8) is not a power series in n. 

The statement in (5.8) is checked to be an identity between formal series (as in the corresponding 
check in [GG2]). 

As a first step to bypass the formal level, the series (5.7) defining M^(x;n) has to be shown to 
be really convergent. This will be true because in the evaluation of M^(x;rf) the only involved 
propagators have scales [p] with p < n — 1 so that, see the factors ip n (x 2 ) , \n (x 2 ) in (5.5), their 
denominators not only do not vanish but have controlled sizes that can be bounded below propor- 
tionally to x 2 by (4.9), i.e. simply by a constant times C 2 \u\~ 2t<) . Using this fact one can actually 
show that the matrices M^(x; rj) are well defined and satisfy symmetry properties similar to (5.4). 

Furthermore cancellations similar to the maximal KAM tori cancellations hold even in this case 
so that MW(i;))) has a special structure, as described in the following Lemma (see Appendices 
A5 and A6 for a proof). 

Lemma 3. Let e < £0 with £0 small enough, and e £ 1(e) — (e/4,e]. Define also rj through 
£ = rj ko . If 1 < n < no, with no defined in (4.. 10)), the following properties hold, 
(i) The series defining the matrices M\- n ^(x; rj), x = u> ■ v, converge and the matrices satisfy the 
same symmetry properties noted for 

EM [ - n] (x;n)E = [M [ - n] (-x; 77)] T , [M [ - n] (x; 77)] * = M [ - n] (-x; 77), (5.10) 
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where E is the 2N x 2N symplectic matrix, (1-4)- Hence, by (5.10), (ix — M\- n ^)E is Hermitian. 
(ii) The matrix M)- n \x; rj) can be written in the form 



where, if Xn-i(x 2 ) ^ 0, 



p[<»](V „\-( d AA H + ^Aa(^ V) d 2 AB H + ^Afl(^; V) \ _ fp[<n]/_. _Nlf 

" ( ' " } " 1 3| A ^ + eM^ (x; n) d% B H + eM§$ (x; n) ) ~ ^ " ^ ^ ' 



Q[ <n ]{x . £)= fixe^2\x; V ) ixeH^ ] (x; V )\ (512) 

/•■•I rri; ■ Tt[ — n ] —[n] . . — [< n] , . — [<n- 1] . . 

(in) lhe entries M i are such that the corrections M 7 1 i{x; r q) = M , \x;r]) — M^^, [x;n) 
are bounded, uniformly in x and e for e small enough, by B e~ Kl2 1 for n < no and for suitable 
no- independent constants B,K\,t > 0; one can take t = t 
(iv) One has 

\\d x M^- n \x, ry)!! < Be 2 , \\d £ [M^ n] (x,r 1 )-M [a] (x,7 1 )]\\<Be, (5.13) 

where the derivatives must be interpreted in the sense of Whitney and the constants B,k\,t > 
can be taken the same as in item (Hi). 

Remarks. (1) The symmetry property (5.10) is proved in Appendix A5. Note that at the first step 
M\°\x;ri) satisfies it, see (4.6). 

(2) The key property in (5.12) is that some entries of Q^- n \x;n) and R\- n \x; rj) are propor- 
tional to x or x 2 : this is proved by exploiting cancellations among families of self energy clusters, 
as described in detail in Appendix A6; note that the single self-energy clusters contributing to 
M\- n ^{x; if) do not have in general the structure in (5.12) and only their sum has. 

A crucial technical point in the proof of Lemma 3 is the fact that, if the scale of a line t is 
smaller than no, as defined in (4.10), then the corresponding propagator admits a bound that is 
qualitatively the same as (4.9). More precisely the following result holds. 

Lemma 4. Let e G (e/4,e]. The propagator g^(x;rj) admits the same bound as g^-°^(x;ri) in 
(4-9). For 1 < n < n , with n given by (4-10), the propagator on scale [n] can be bounded by 

n] (x;v)\<^, (5-14) 
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for some positive constant C. 

The proof of Lemma 4 proceeds as that of Lemma 2, in Appendix A4, and we do not repeat 
it here. Once the bound (5.14) is established, the proof of convergence is the same as the one 
discussed in Appendix A3 of [GG2] and we do not repeat it here. Item (ii) of Lemma 3 simply 
follows from convergence and from the remark that M\- n \x; n) — Ai^(x;rj) is of order e 2 (cf. 
Remark (3) after Definition 4). The bounds in items (iii) and (iv) of Lemma 3 also follow from 
the proof of convergence, see Appendix A3 of [GG2]. 

We have therefore constructed a new representation of the formal series for the parametric 
equations for the invariant torus: in it only trees with lines carrying a scale label [— 1], [0], . . . , [n — 1] 
or [> no] and no self-energy clusters are present (note that, so far, self-energy clusters may have 
only scales [n] with n < n ). The above lemma will be the starting block of the construction that 
follows. 



6. Renormalization: the infrared resummation 

From the proof of Lemma 2 it is clear that for x 2 < pn 2k °^ 1 and n > n it will not be possible to 
bound g\ n \x; n) by a constant times x~ 2 . 

So, the first problem to face when reaching scales n > n is the computation of the eigenvalues 
of (ix — M\- n \x;r]))E, in terms of which an estimate of the size of the propagator g\ n ^ can be 
deduced. 

If M\- n \x;rf) does have the structure in (5.12) then an approximate computation of the eigen- 
values of (ix — A4^- n ^)E leads to the following Lemma, proved in Appendix A7. 

Lemma 5. Let n > n and let us assign a matrix M\- n ^(x;-q) satisfying the properties in (5.12) 
admitting right and left derivatives with respect to x and e, bounded as the derivatives in (5.13), 
and having the structure described by (5.12), with the entries M^™}(x;r)) such that the correc- 
tions m!™L(o;; rj) = M^3(x\rj) — m]^Z ^(x^rf) can be bounded (as in item (iii) of Lemma 3) 

by ,(x\rj)\ < £>e~ Kl2 " /T1 , for suitable constants B, m and T\. Then the uniform norm of 

(ix — M\- n }(x; 7])) can be bounded below by 

\\ix-M^(x;r))\\ > ^- min{x 2 ,\x 2 - X^(x; V )\}, (6.1) 

where fi = jJ^ is the largest eigenvalue of and X^(x;rj) — l^(x;n)eri ka ^ 1 is a real function 
with £^(x;r)) — ck d 2 3 H a f3 , [ a _1 (l + 0(rj)). Furthermore X^(x;n) is right and left differ entiable 
in e and x and the derivatives satisfy the following dimensional bounds: 

C-i^o-i < \d±\W(x;r])\ < Cn k °-\ |^A w (x; r])\ < Cn k ° , (6.2) 

for some positive constant C. 

Remarks. (1) The bound (6.1) suggests to replace the classical small divisor x 2 used in previous 
sections by the (n-dependent) quantity min{x 2 , \x 2 — X^ n \x; rj)\}, that is essentially what we shall 
do below. In particular we shall replace the argument x 2 of the support functions ib, x by a quantity 
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A^(x; n), behaving like minjcc 2 , \x 2 — X^ n \x; which will be the measure of the strength of the 
resonance for scales larger than n . With this choice we shall in general introduce a singularity 
in the definition of the propagators and self-energy matrices on scales n > uq: this is due to the 
presence of a minimum and of an absolute value in the definition of A^(x; rj). This is ultimately 
why in (6.2) the right and left derivatives appear, rather than the plain derivatives, as in (5.13) 
above. This could be avoided by using a smoothed version of the quantity A^{x;rj) introduced 
below, but we shall not discuss it here. 

(2) The bound on the derivatives of X^{x;rj) with respect to e follows from the expression of 
M\^{rj), in (4.3) to (4.5), and from the bounds (5.13), which allow to control the corrections. On 
the contrary the bound on on the derivatives of \^(x; rj) with respect to x does not follow directly 
from (5.13), and it is explained in Appendix A7. 

Depending on the sign of e, from now on, the analysis changes qualitatively. If e is such that 
cePi"^ 1 < (so that X^ n \x; rj) in (6.1) is negative), the minimum in (6.1) is always realized by the 
classical small divisor x 2 . This implies that in this range of scales it is possible to proceed in the 
same way discussed above in Section 5. We do not repeat the details, and we concentrate on the 
opposite case, namely e with the sign such that X^(x; rj) > 0, which presents new difficulties. 

In this case convergence problems can still arise from the propagators gi- n °\ (x; 77), which become 
uncontrollably large for x = u> ■ v close to the eigenvalue X^ n \x;rj) in (6.1). 

We introduce a sequence of self-energies A'™' (n) , [GG2] , representing the locations of the singu- 
larities of ix — M\- n \x;vj) and, correspondingly, we introduce a sequence of propagator divisors 
At"! {x; rj), which will give a bound for the size of the propagator on scale n > no. Then, we modify 
the scale decomposition, measuring the strength of the singularity in terms of A^ n \x;rj) rather 
than, as done in Section 5, in terms of the classical small divisor x 2 . 

The choice of the scale decomposition will be done in such a way that the dimensional bounds 
for the propagators will be the same as for those with scales n < n : this will reduce the analysis 
of the infrared resummations to the same convergence proof discussed in Appendix A3 of [GG2] . 

We introduce the following definition. 

Definition 7 (Self-energies and propagator divisors). Let the function X^(x;n) be as in Lemma 
5. 

(i) The sequence of self-energies }^ n \rj) is defined for n > uq by 

A N („) d i! AN (^/a [ ™- 1] (r)),v) , A'™ ' (rj) d = A'™ ' (0; rf), (6.3) 
provided XS n \r]) > 0, n > no. 

(ii) The propagator divisors are defined for n > no by 

A^{x;n) d = f min{x 2 , \x 2 - X [n] (n)\}. (6.4) 

By repeating the analysis of Section 5 we can represent the function X(tp) via sums of values 
of trees whose lines can carry scale labels [— 1], [0], . . . , [no — 1], [no], [no + 1], . . . and which contain 
no self-energy clusters (i.e. they are renormalized trees; see Definition 5 in Section 5). The new 
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propagators will be defined by the same procedure used to eliminate the self-energy clusters on 
scales [n] with n < n — 1. However for n > n the scale of a line will be determined in terms of 
the recursively defined A^{x;rj), (6.4), rather than in terms of x 2 , see (5.2). 

Set X^Cr) =' YC=oXm(x 2 ), Y n {x;ri) d ^ Ul=n Xm{* [m] {x;ri)) for n > n and F„ _i = 1. 
The definition of the new propagators will be 

ff [n„ + l] "e/ Xno _ i(a;) (A [n ] {x . v)) ^ ng + l{A [no + l] {x; ??)) (fl .2 _ ^[<n + l] (x; ^ ^ ^ 

and so on, using indefinitely the identity 1 = ip n (A^(x; r])) + Xn{A^ n \x; rj)) to generate the new 
propagators. 

In this way we obtain a formal representation of X(tp) as a sum of tree values in which only 
rcnormalized trees appear and in which each line £ carries a scale label [ne]. This means that we 
can formally write X(ip) as in (5.8), with Val(0) defined according to (5.9), but now the scale label 
[ne] is such that ng can assume all integer values > — 1, and no line carries a scale label like [> n]: 
only scale labels [n] are possible. 

We can summarize the discussion above in the following definition. 

Definition 8 (Propagators and self-energy matrices). Given a sequence of 2N x 2N matrices 
M^ m \x;r]), m > 1, let M^{x;r]) = M^ n \x;r)) - M^- n ^(x; n) with M^(x;r)) = M^(n) 
(see (4.2)), so that M [ ^ n] {x;n) = ^„M H {x;r}). Setting A^(x;t]) = x 2 if n < n , define for 
all n > 

M^ ln] (x; V ))U n m >oXm(A^(x: V )) 
9 M = x 2 -M^(x; V ) • (6 ' 6) 

(for n — this means ipo(x 2 ) {ix 2 — M^{ri))~ l ). We say that = g^- n \u> ■ Vf,e) is a propagator 
on scale [n]. The matrices M\ m \x; rj) will be defined as in Section 5 for n < no and will be defined 
recursively also for n > no in terms of the self-energy clusters S^ n l introduced in Definition 4, 
Section 5, setting for n > no (cf. ( 5. 7) ) 

M^(x;n) = (f[ Xm (^ m \x;f 1 )))f2 £ V T (x;n), (6.7) 

m=0 k=2 T£S]^ ri _ 1 

where the self-energy values VT(x;n) are evaluated by means of propagators on scales less than [n]. 
Note that we have already defined (consistently with (6.7))) the matrices M)- n ^ with n < n and 
the propagators on scale [—1], [0], . . . , [no — 1] (so that (6.6) defines also g^ n °\x;rj)). 

Of course the above definition makes sense only if the series (6.7) can be shown to be convergent 
for all n. For this purpose an inductive assumption on the propagators on the scales [m] , < m < n 
is necessary. 

Inductive assumption. Let no be fixed as in Lemma 3. 

(i) For < m < n — 1 the matrices M.^ m \x;r]) are defined by convergent series for all \e\ G /(e) = 
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[e/4, e] and, for all x, they satisfy the properties (i) and (iv) of Lemma 3. Moreover they can be 
represented as in (5.11), with p[-"l , Q\S n \ , i?[- n l as in (5.12) and the entries My"' bounded as 
described after (5.12). 

(ii) There exist K > and open sets £° l7 m = 0, ...,n, with £° n C /(e), smc/i £/ia£, defining 
recursively A'"^(n) in terms of A' m_1 ' {if) for m = no, . . . ,n — 1 fey (%) m Definition 7 above, while 

setting A' m ' (n) = for m = 0, . . . , no — 1, and defining n d = r (l + 5i) + iV, wit/i <5i > 0, one /ias 
/or e £° 

rH(, ; „)^ 1x1-^(77) — im c ° 



> 2"5' 



f T1 



with e = rj ka , for all m < n — 1 and a?/ x. 



(6.8) 



Remark. As in [GG2], a key point is checking that A^(x;n) can be used to bound below the 
denominators of the non- vanishing propagators on scale [n]. If x has scale [n], with n > n , one 
has 

x 2 - AM (a:; n) > x 2 - A w (77) - A [n] (n) - A [n] (x; r?) 



1 

> - 

- 2 

1 

> - 

- 2 



•A [nI (r/) +2-(™+ 2 )C - X [n] {-\f^~%),v) ~X [n] (x;v) 
A [n] (v)\ \\x 2 -M^(x, V )\\>-^AM(x;r,), 



(6.9) 



having used the lower cut-off t/} n (A^ (x; rf)) in the propagator (sec (6.5)) to obtain the first two 
terms in the second line while the upper cut-off Xn-i(A^- n ~ 1 } (x; 77)) has been used to obtain posi- 
tivity of the difference between the second and third terms in the second line, after applying the 
second inequality in (6.2), so that the last term in the second line of (6.9) can be bounded above 
proportionally to e2~™Co. 

The inequality (6.9) allows us a complete word by word reduction of the proof of the inductive 
assumption above to the corresponding inductive assumption of [GG2]. The symbols here and in 
[GG2] have been chosen to coincide so that the analysis in Section 6 and Appendix A3 of [GG2] 
can be taken over and reinterpreted, with no change, as proofs of the above inductive hypothesis, 
apart for the check of the measure £^ in (6.8), which in the present case is slightly different from 
the corresponding computation in [GG2]. The estimate (6.8) of the measure of £^ is explicitly 
given in Appendix A8. 

We can summarize the previous discussion into the following Lemma. 

Lemma 6. There is en small enough such that for e < £0 and e G /(e), if the inductive hypothesis 
is assumed for < m < n — 1 then it holds for m = n. 

The series for X(ip) is now fully "renormalized" and its terms are well defined for e in a set £ 
whose measure is large (near 0). The series is expressed as a sum of renormalized tree values of 
tree graphs without any self-energy graph. Therefore the series is convergent (by Sicgel's lemma); 
see the corresponding discussion in Appendix A3 of [GG2]. 
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In the derivation geometric series with ratio z > 1 have been considered and the rule Y^k=o z>t = 
(1 — z) _1 has been repeatedly used: this is not mathematically rigorous and therefore an a posteriori 
check must be made that the function X(ip), via (1.5), actually does satisfy the equations (1.4). 
The check, however, is a repetition of the corresponding (essentially algebraic) check in Appendix 
A5 of [GG2]. 

Therefore the proof of Theorem 1 is complete 

7. Conclusions. 

We conclude by mentioning some problems that seem to us of interest. 

(1) First we note that the proof of Theorem 1 also provided some informations about the analyt- 
icity properties in 77, hence in e, of the surviving lower-dimensional tori; cf. [GG1] and [GG2] for 
further details. As in the quoted papers, it can be interesting to further investigate such properties, 
and related ones, such as Borel summability. 

(2) The uniqueness of the resonant tori appears to be a hard problem. Even a proof of Borel summa- 
bility would not resolve the problem as there could be solutions which are not Borel summable. 
Just as in the maximal tori case analyticity is not sufficient to guarantee uniqueness. The very 
recent [BT] does not settle the question, as it proves uniqueness of the C°° diffeomorphism map- 
ping the unperturbed invariant tori which are conserved into those which are explictly constructed, 
but it dose not eliminate the possibility that other nearby quasi-periodic motions with the same 
rotation vectors exist. 

(3) Another open problem concerns the possibility of removing the assumptions we made in this 
paper. We have been able to do this in several particular cases but we did not find a satisfactory 
general formulation: for instance can something be said in general in the case in which the average 
of the perturbation vanishes identically? 

(4) Finally, it would be interesting to understand what happens in the case of n-dimensional tori, 
with n strictly less than N — 1 (that is in the case of several normal frequencies), by relaxing the 
assumptions on the perturbing potential with respect the results existing in literature, as done 
here and in Cheng's paper [Chi] and [Ch2] for n = N — 1. This is a substantially more difficult 
endeavour with respect to that considered here and in Cheng's quoted papers. 

Appendix Al. Heuristic analysis leading to fractional series 

Consider the simple case H(A 7 B,a,0) — oj ■ A + 5 (A 2 + B 2 ) + ef(a, 0), and let assumptions (a) 
and (c) in Section 1 be satisfied by /. Let us consider the generating function 

T(A', B', a, (3) = a ■ A' + 0B' + e*i(A', B', a, 0) + e 2 V 2 (A\ B', a, 0)+s 2 a-C + e 2 0p, (Al.l) 
where, if A d = a; • d a , 

*i(A',B',a,/?) = f -A- 1 [l-(A'-9 a + B'd ) A" 1 ]/^, 0), (7.2) 
and, if $(A', B', a,p) = \ [(c^) 2 + (^j*i) 2 ] , 

* 2 (A',S',a,/3) = f -A- 1 [l-(A'-9 a + B'dp) A- 1 ]$ #0 (A', B', a, 0), 
C = -0A'*o(A',O,A,)| A , =o , p=-d B «f>o(O,B',0 o )\ B , =o . 
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The canonical map (A,B,a,0)< — >(A', B',a',0') generated by T(A' , B' ,a, 0) transforms the 
Hamiltonian H(A, B, a, 0) into 

H'(A', B', a', 0') = w • A' + i(A' 2 + B' 2 ) + sf o (0')- (A1A) 
- e 2 dp,f (P') ■ A- 2 ^/ /0 K, 0') + e 2 $o(0, 0, 0') + 0(e/' 2 ) + 0(e 3 ) . 

Writing the Hamiltonian equations generated by (Al.4), we realize that in a small neighborhood 
of /3o the evolution equation for 0' takes the form: 

0' = -ce(0' - a ) k ° - e 2 a + e 2 a\0 - O ) + O{s(0' - o ) ka+1 ) + O(e 2 (0' - 0o) 2 ) + 0(e 3 ) , (A1.5) 

where a = dpi$o(0, 0, 0o) is the same constant defined in (1.9) and a' is a suitable constant. 

Therefore, under the assumption that both a and c are non-zero, the perturbed equilibrium of 
the angle 0' is equal, up to high order corrections, to 0' Q = 0o + S0o, with 

( (-ae/c) 1 /^ if fc is odd, ( . , 

' " \ ±(-ae/c) 1 /fco if fc i s even and ae/c < 0. k ' 

If fco is odd and ce > the approximate perturbed equilibrium point 0' o is quadratically stable, 
while if ce < 0' o is quadratically unstable: hence we shall say that the resonant invariant torus 
is elliptic if ce > and hyperbolic if ce < 0. The center of oscillations is displaced by 0{e x / ka ) so 
that a fractional series in powers of r] = e 1 ^ has to be expected (at best). 

If fc is even the unperturbed equilibrium point O can be continued into a perturbed equilibrium 
0'o = 0o + $00 (with 50o vanishing as e — > 0) only if aec < 0. In this case, if ce50o > the 
approximate perturbed equilibrium point 0' o is quadratically stable, while if ceS0o < 0' o it is 
quadratically unstable: hence we shall say that the resonant invariant torus is elliptic if ce50o > 
and hyperbolic if ceS0o < 0. Given the results known for elliptic and hyperbolic resonances 
Theorem 1 becomes a natural conjecture. 

If c 7^ and a = the theory depends on the corrections in (A1.5): for instance if a' ^ a 
natural conjecture is that the new equilibrium point for 0' is approximately 0' o + S0o, with S0o 
satisfying the equation (S0o) k °~ 1 — ea'/c, whenever this equation is solvable. The stability of 
this equilibrium point can be again analyzed in terms of the relative signs of e, a, c. If a' = one 
expects that the higher order terms have to be studied in details. 

Appendix A2. On assumption (a) 

Here we want to show that assumption (a) is not a purely technical assumption; that is we want 
to shows that generically a Hamiltonian of the form (1.3) violating assumption (a) cannot admit 
quasi-periodic motions of codimension 1 continuously connected to an unperturbed motion of the 
form (1.6) in the limit £ — ► 0. We show this by producing an explicit example, given by an 
Hamiltonian of the form (1.3), satisfying assumptions (c) and (b) and violating assumption (a), for 
which we can prove absence of quasi-periodic motions of codimension 1 tending to an unperturbed 
motion of the form (1.6) as e —> 0. 
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The counterexample is given by the following Hamiltonian 

H = A+\A 2 + \b 2 +e ( / ^^+sin/3cosa) , (A2.1) 

where A, B £ M 1 and a, (3 G T 1 . If e = the Hamiltonian (A2.1) is strictly convex (then it satisfies 
assumption (b)) and it admits the unperturbed periodic motions 

A{t) = 0, a(t) = 1, B{t) = 0, 0(t) = A), (-42.2) 

parametrized by the choice of (3q G T 1 . Choosing /3 = 0, we see that /o(0) = 9/3 /o(0) = 9|/o(0) = 
and 9^/o(0) = 2 ^ (so that assumption (c) is satisfied with k = 2). Using the fact that / 
is independent of the action variables and that f(a, /3 = 0) = dpf(a, (i = 0) = 0, we see that the 
quantity a defined in (1.8) is identically 0: this means that assumption (a) is violated in the case 
under analysis. 

We now investigate the possible existence of periodic motions of the form 

a(t) =t + a(t;e), A(t) = a(t; e) , (3{t) = 77(e) + b(t; e), B(t) = b(t; e), (A2.3) 

where 77(e) is a continuous function of e such that lim e ^ v( £ ) = and a(t; s),b(t; e) are average 
periodic functions of t of period 2ir. We want to show that it is impossible that a motion of the 
form (A2.3) satisfies the Hamiltonian equations of motion 



(A2A) 



(A2.5) 



d(t; e) = e sin(i + a(t; e)) sin(?7 + b(t; e)) , 

b(t; e) = — e[sin 2 (?7 + b(t; e)) cos(?7 + b(t\ e)) + cos(t + a(t; e)) cos(?7 + b(t; e))] . 
In fact a average periodic solution to (A2.4) is necessarily of the form 

a(t;e) = — e(r/sint + ^ sintcost) + e 2 0(e 2 + if), 
b(t;e) = ecost + eO(e 2 + rf). 

Averaging over t the second of (A2.4) we find 

/ — (cos(i + a(t; e)) cos(?7 + b(t; e)) + sin 2 (r; + 6(t; e)) cos(?7 + b(t; e))) = 0. (A2.6) 

JO 27T 

Now, using the expressions (A2.5), we see that (A2.6) is equal to ?7 2 + plus terms of order at 
least i] 4 + e 4 , and this leads to a contradiction. 



Appendix A3. Counting the number of trees lines 

We want to show that the number of lines of a tree contributing to A^ k \ B^ k \ oS k ^ or b^ can 
be bounded, above and below, by an O(l) constant times k. A lower bound k/ko is an immediate 
consequence of the definitions. 
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We proceed by induction, using that: 

(1) the trees contributing to A( fc °) and B^ k °^ have one line (see (2.2)); those contributing to A( feo+1 ) 
and 5( fco + 1 ) have two lines; 

(2) the trees contributing to and b^°\ with v ^ 0, have at most two lines and those con- 
tributing to a!£° +1 ^ and b^ a+1 \ with u ^ 0, have at most three lines; 

(3) = Pi is represented by a trivial tree (i.e. a tree with one line) and the trees contributing 

(2) 

to b ' = fa have at most ko + 2 lines. 

We introduce the following definitions: 
Bk{I) is the maximum number of lines of the trees contributing to A^ k \ B^; 
-Bfe((/?) is the maximum number of lines of the trees contributing to b^\ u =/= 0; 
Bk(b ) is the maximum number of lines of the trees contributing to b^ = (3k- 

We shall make the following inductive assumption: 



B k (I) < 3fc (fc - fc + 1) - 4fc - 2 

Bk(<f) <3fco(fc-fc + l)-4fc -l 
B k _ ko+1 (b ) < 3fc (fc - k + 1) - 4fc 



k > k Q + 1, 
k > k + 1, 
k > k + 1. 



(A3.1) 



By the remarks above and the fact that fc > 2 it follows that at the first step (that is k = fc + 1) 
the inequalities above are verified. 

Assume inductively the inequalities in (A3.1) for fco + l<fc</i — 1 and let us prove them for 
k = h > ko + 2. Let us start with the third inequality. We call Co the first node preceding the root. 
By the rules explained in Section 3 we have that the sum of the orders of the s = s Vo subtrees 
entering Do must be h. So, using the inductive assumptions, and calling (see Figure 3): 



si the number of subtrees of type I^ k \ k > ko + 1, entering Co; 

s'j the number of subtrees of type 1^ entering Co; 

s v the number of subtrees of type tp v , k > ko + 1, v ^ 0, entering Co; 

s' v the number of subtrees of type (f~v°\ v ^ 0, entering Co; 

s the number of subtrees of type b Q k \ k > 2, entering Co; 

s' the number of subtrees of type b^ entering Co, 
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si 
s'i 

<p 

so 
s'o 

Figure 3. The lines entering i>o represent symbolically the bundle of subtrees en- 
tering Co and with root lines of types I^ k \ I^ k °\ ip^\ ip^ a \ with k > fcrj + 1, and 
/3fc, /3i, with k > 2. The number of subtrees in each bundle is indicated by the right 
labels. 

Bh-k +i(bo) can be bounded by the maximum over the choices of sj, s'j, s v , s' v , s , s' of: 

l+3A;o/i-{(3^+A;o+2)s / +(3^-l)s^+(3^+A;o+l)s v +(3fc^-2)s;+4fcoSo+(3A;o-l)s }. (A3.2) 
Let us first consider the case s/ + s'j + s v + s' > 1. 

In this case, if s = 1, since /i > fco + 2, we must have that either s v or sj is = 1, and, in both cases, 
we can bound (A3. 2) by 1 + 3koh — (3&Q + k + 1), that is the desired bound; if, on the contrary, 
s > 2, then (A3. 2) can be bounded by 1 + 3k h - (3fcg - 2) - (3fc - 1) < 3k h - 3fcg - fc . 
Let us now consider the case si + s'j + s v + s' = 0, in which case s = so + s' Q and necessarily s > ko 
(note that in this case the s + 1 derivatives in (3.6) must be derivatives with respect to (3 so that 
s > ko). If s > fco + 1, then (A3. 2) can be bounded by 1 + 3k h— (fc + l)(3fc — 1) < 3fc /i — 3fcp — fc . 
If s = fco, by the rules in Section 3 it must be s' < ko — 2, so that (A3. 2) can be bounded by 
1 + 3k h - 4fc (fc - sd) - (3fc - l)s < 1 + 3fc /i - 4fcg + (fc - 2)(fc + 1) that implies the desired 
inductive bound. 

Let us now consider the first of the three inequalities in (A3.1). By the rules explained in 
Section 3 we have that the sum of the degrees of the s subtrees entering Do is h — ko or h if, 
respectively, <5o = 1,0. If <5 0o = 1, Bh(I) can be bounded by the maximum over the choices of 
Si, s'j, s v , s' v , so, s of 1 + 3ko(h — ko) — {■ ■ ■}, where the brackets include the same expression as 
in (A3. 2), and we can bound this number by 1 + 3koh — 3fcg — (3fc — 1) that is even better than 
the desired bound. If <5 0o = 0, one must have s = s/ + s'j > 2 (as s = 1 corresponds to a not 
allowed tree graph) and Bh(I) can be bounded by 1 + 3fco^ — 2(3/cq — 1), that is even better than 
the desired bound. 

We finally consider the second of the three inequalities in (A3.1). By the rules explained in 
Section 3 we have that the sum of the orders of the s subtrees entering Do must be h — ko or 
h, if, respectively, <5 0o = 1,0.. If S Vo = 1 the number of lines of the corresponding trees can be 
bounded in the same way as we proceeded above for the first inequality in (A3.1). If <5 0o = 0, only 
lines of type / can enter Do, and then the corresponding trees have a number of lines bounded by 
1 + 3k h - {(3kl + ko + 2)sj + (3fcg - l)s'j}. So, if s/ > 1 the desired bound follows. If s/ = 0, 
since h > ko + 2, we must have s'j > 2, and again the bound follows. 
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Appendix A4. Proof of Lemma 2 

To estimate the eigenvalues of (ix — A4^ (ri))E, we start by rewriting (ix — (r\))E in the form 

(i X -M^))E=( ix - + R Q , - p +Q ). (A4.1)) 

Note that, by the strict convexity of the Hessian of Ho if e is small enough, P admits N positive 
eigenvalues < /4°^ < • • < fj,jp = fx. Let v a (x;rf) = (^^j be an eigenvector of (ix — M^)E 
with eigenvalue a (here u and v are two column vectors of dimension N): 

-R -ix + Q\ ( u\ _ ( -Ru + (-ix + Q)v\ _ (u 



ix + Q^ P )\v)~\ (ix + Ql)u + Pv ) a \vJ- (A4 - 2) 

Now, either \a\ > Mi^/^ or -j < P — a < In the latter case (P — a) is invertible, 2/(3/^) < 
(P — a) -1 < 2/[if \ and we can rewrite (A4.2) as 

v = -(P - a)- 1 (ix + Q^)u, -Ru-(-ix + Q)(P-a)~ 1 (ix + Q t )u = au. (A4.3) 

The second equation in (A4.3) is of the form 

x 2 (P - a)- 1 + O(ex) + 0(erj k °- 1 )] u = au, (AAA) 

so that | a | > j^x 2 + 0(ex) + O(er] ko ~ 1 ). If x 2 > p\e\f] k °^ 1 , for p large enough, the latter estimate 
implies \a\ > j^x 2 , and Lemma 2 is proven. 

Appendix A5. Symmetry properties of the self-energy matrices 

In this section we discuss the symmetry properties (5.10) of M\- n \x; rj). 

We begin with proving (5.10) for M\ n \x;rj). We inductively suppose that the same symmetry 
properties hold for both Ai^(x;r]) and g^(x;rj), for < p < n, 

Egb\ X ;rj)E = [g [p] (-x:ri)] T , ]gW ( x - v )]* = ff M (-x; rj), (Ab.l) 

where the * denotes complex conjugation and T transposition. Note that if p = 0, M.^(rj) and 
(x; rj) satisfy the desired symmetry properties, as discussed in Section 4. 
Consider the representation of M\ n \x]ij) given by (5.6) and (5.7): in the following confusion 
between T denoting transposition and T denoting a self-energy cluster will be avoided by renaming 
the cluster T with a new symbol D. Note that F*(u$) = F(— u$), so that, using the second of 
(A5.1), [Vd(x;ti)]* can be written as 2 

[VD{x;v)r = ( ^wmr( n eS °)( n f *(»»)){ n i^h^vw) = 

' ^ veV(D) vev(D) leA(D) 

m) WD *r jj £ s„)( jj F{ _^( jj g M { - Xt .,r,))= VD .{-x; e ), 



\MD)\ 



We stress that we use the convention that the internal 7 indices of D are summed over, while the momentum 
labels are not summed over; in this way the value Vrj> explicitly depends only on: the momenta of the internal nodes, 
the external momentum and the two 7 labels associated to the incoming and outcoming lines. 
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where D* is a cluster topologically equivalent to D, with opposite mode labels (clearly the corre- 
spondence D< — >D* is one-to-one). Summing over the choices of D (see (5.7)) the second of (5.10) 
follows. 

To prove the first of (5.10), it is convenient to shorten notations; see Figure 4 for a pictorial 
representation of the symbols. Given a self-energy cluster D, let us denote by V and V' the nodes 
such that the exiting line of D comes out from V' and the entering line of D enters V. 

(1) Let V' = Vi, . . . ,V n = V be the nodes on the path C = C(D) joining V to V' , and Vj be their 
mode labels; let also lj = (Vj+iVj) be the line joining the two successive nodes Vj+i and Vj in C; 
let 7^ , 7j be the component labels at the beginning of the line in C exiting from Vj and at the end 
of the line in C entering the node Vj (they are, respectively, the lines lj-i and lj). Recall that, by 
definition of self-energy cluster, any propagator gj associated to the lines tj G C has a scale [n<J, 
with ni j > 0, i.e. there cannot be 0-momentum propagators along the path C. 

(2) Let -dj be the (possibly empty) family of subtrees of D with root in Vj and with no line in 
common with C (not represented in Figure 4: i9j can be imagined to be a set of trees with the 
root line ending in Vj or, possibly, ftj may consist just in Vj). 

(3) Let Nj be the N x N symmetric matrix dy^.f^ (/3o,#j), where f Vj ((3o,'dj) is a function 
depending only on the set of trees i?j and on f v (0) (it can be read off (5.6) and (5.7)) and <9y, <9 7 
must be interpreted as explained after (3.2). 

(4) The momentum flowing in a generic line I of the graph D will be the sum of the momentum 
that would flow on the line if the entering line momentum v was plus v if the line is on the 

path £; then the propagator matrix gj associated to the line lj of £ is equal to g^ le i\xj + x;rj), 
where is the scalar product of u; • and x = u> • v. 



Figure 4. A self-energy cluster D and the path C connecting its external lines. The 
subtrees internal to D with root on C arc not drawn. 



Then the matrix Vd(x; rj), see (A5.1), can be concisely written as 

(Vz>(a;;»7)) 7 , 7 = (EN m EN 2 g 2 ...g n _ 2 EN n _ 1 g n _ 1 EN n ) . (A5.3) 

Interchanging the external lines (i.e. having V 1 as entering node and V as exiting node) generates 
a new self-energy cluster D 1 in which the momenta flowing in the lines of the subtrees ftj with 
roots on the nodes Vj G C are unchanged while the momentum on the line lj £ C changes from 
vf. + 1/, with equal to the momentum which would flow on lj if the external momentum v was 
set equal to 0, to a new value — 1/° + v. 
The matrix Vre(x; rj) can therefore be written 

(V D >{x; V ))^ = (EN n g' n _ 1 EN n _ 1 g' n _ 2 ...g 2 EN 2 g' 1 EN 1 ) i , 7 , (A5A) 

where g'j = g^ n *i\—x®, + x;rf). Inserting 1 = E T E after N n ,...,N 2 and using the symmetry 
of Ni, . . . ,N n , the inductive validity of the first of (A5.1) and the transposition rules for matrix 
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products, we get 

(V£,(s; v)) Yl - (NrflEN 2 & . . . g'^_ 2 EN n ^ig'^_ l EN n E T ^j , (A5.5) 

where g" = g^ nej \x®. — x;rf). Hence EVp,(x;ri)E = Vd{— x;rf) completing the inductive proof of 
(A5.1), because M\ n \x; 77) is defined as a sum over D of the values Vd (see (5.7)). 

Remark. From (5.10) and (A5.1) follows that gl- n \x;r]) satisfies the same symmetry properties as 
MW(x;r)): 

Eg^(x;r ] )E=[g^(-x; V )] T , [g^(x; r?)]* = g^ n \-x; r,). (45.6) 



Appendix A6. Cancellations 

In this section we discuss the cancellations needed to show that M\- n \x\rj) has the structure 
described in (5.11) and (5.12). In particular we want to show that the elements (A4^ (x; ffj) , 
with either 7' = Ai or 7 = aj are proportional to x, while the elements with 7' = Ai and 7 = etj 
are proportional to x 2 , for x 2 < prf^^ 1 . 
As in [GG2], the contributions to M\ n \x;ri) coming from clusters such that 

£ Ki>(c /2 6 M) 1/r ° (A6.1) 

OGV(T) 

require no cancellations: in fact the exponential decay as u v -^00 of the node functions implies that 
the contributions coming from clusters satisfying (A6.1) are smaller than Cx 2 , for some constant 
C > 0. 

When (A6.1) does not hold we need to exploit cancellations, because a priori there is no reason 
for some of the entries of M\ n ^ (x; rj) to be proportional to x or x 2 for small x. Following the 
same strategy used in Appendix A2 of [GG2] (and of Appendix A4 of [GG1]), the necessary 
cancellations can be checked to occurr by collecting clusters violating condition (A6.1) into families. 
Such cancellations are due to the same mechanism pointed out first in [Gl], [G2] and [CF]. The 
following analysis follows [GM]. 

Given a self-energy cluster D on scale [n — 1] , let us call D the connected subset of D containing 
no line on scale [—1] and containing the nodes and 0' to which the entering and exiting lines 
are attached. By definition ^2 veDg fx> = 0- Note also that, if the path C connecting and t>' 
is non-trivial (i.e. 7^ t>'), C is completely contained into Do; moreover, if D does not contain 
lines on scale [— 1], then Dq = D. We define the family T = To as the set of self-energy clusters 
obtained from D by shifting the entering and exiting lines by reattaching them to the nodes of Do 
in all possible ways. We then consider the sum 

VAx;v) = £ Vd(x; V ), 046.2) 

contributing to the r.h.s. of (5.7). When summing Vr(x\rj) over all distinct families of clusters 
on scale [n — 1], we recover the quantity M^ n \x; rj) defined in (5.7). To prove that A^-^ has the 
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structure described by (5.11) and (5.12), we proceed in the following way. We first note that, by 
the same analysis of Section 6 of [GG2], Vyr(x;n) is differentiable with respect to x in the sense of 
Whitney; then, we show that (Vy(0; Tf)) , = if 7' = Ai or 7 = ctj and (d x V^(0; T])) , = if 
7' = Ai and 7 = ctj. 

Let us begin with proving that (Vy-(0; 77)) , = if 7' = ^ or 7 = ex,-, and to be definite let us 
assume that 7 = ctj for some j = 1, . . . , N — 1. Recall that, employing the notation of Appendix 
A5, we can write each matrix Vd(x;t]) appearing in the sum (A6.2) in the form (A5.3). By the 
explicit expression of (A5.3) and of the matrices Nj, we see that, when computing Vd{x;t]) in 
x = 0, we find (Vd(0; r?)) 7 ' Qj = 0{F, t>', 7])y (i^t>)j for some function 0{T , t)', r])y independent of 
the choice of the node to which the entering line with label 7 = ctj is attached to. 

If we change D within T by attaching the entering line to all possible nodes in D (keeping 0' 
fixed), we see that the global contribution from such self-energy clusters is equal to 0{T , 0', r?)y 
• J2 oeDo {vo)j, which is by the property ^2 oeDo (^o)j = 0. So, if 7 = ctj, the proof of the fact 
that (V;f(0; J?)) , = is complete. In the case that 7' = Ai the proof is completely analogous and 
we do not repeat it here. 

Let us now turn to the proof of the fact that, if 7' = Ai and 7 = ctj, then (d x Vj^(0; rfj) , = 0. 
First note that, by the explicit form (A5.3) of the value Vd{x;t]) and in particular by the fact 
that Vd(x;t]) depends on x only through the propagators along the path C, we find that, when 
differentiating Vd(x; rj) with respect to x, the effect of the derivative is as follows: 

n-l 

d x V D (x;r]) = ^ (eN i9i ■ ■ ■ g i - 1 EN i }d x g i (EN i+1 g i+1 ■ ■ ■ g^EN^. (A6.3) 

i=l 

Given the line li E -Do, we call D\ and D 2 the two connected distinct subsets of D n obtained from 
Dq by detaching the line £i and such that D\ contains the node 0' to which the line exiting from 
D is attached, while D2 contains the node to which the line entering D is attached. D\ and D2 
are two subgraphs of D with two external lines, one coinciding with one of the external lines of D, 
the other coinciding with £i, see Figure 5. 




Figure 5. Subdiagrams D\ and D2 in a self-energy cluster D. The drawn lines 
connected to and 0' are the external lines of D. 

If we define the values of £>i,Z?2 by formulae analogous to (A5.3), we can rewrite each term 
under the sum in (A6.3) as V Dl {x;v)d x g ln ^ ] {x° i + x; 

t))Vd 2 {x; tj), where, with the same notations 
of Appendix A5, we defined x® = u> ■ v®. , with 1/°. equal to the momentum which would flow on ii 
if the external momentum u was set equal to 0. 

In particular let us consider the value (Vdi (0; fj)d x g Lnt ^ (x°; t])Vd 2 (0; 7 7))A i a J at x = 0, that is 
one of the contributions we are interested in, and again note that by the very definition of Vdx 
and Vd 2 , it holds that (Vd^O;^)) = (iv)iQ(ro) 7l and (Vd 2 (0;v)) = ^( ro ')72( I> 'o)i, where 
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tT>,tu' are the two nodes where ii enters into and exists from, see Figure 5, and 0(tt)) 7l , f2(tt)') 72 
are two functions depending on the structure of D x and D 2 , but not on the choices oft)', t) within 
D U D 2 . 

This implies that if we change D within T by attaching the entering line to all possible nodes 
in D 2 n D and by attaching the exiting line to all possible nodes in D\ n Do, keeping £ t fixed, 
the sum of the contributions of the form {Vd 1 (0; r])d x g^ ne ^ (x°; ti)Vd 2 (0; ?y)) J 4 i Q j from this class of 
graphs is equal to 

( J2 M^M,Ad x g lne '\^;v)h 1 ,M^% 2 ( E ( A6A ) 

ceflinflo veD 2 nD a 

Note that, given any graph D with the structure in Figure 5, the graph D' in which the entering 
and exiting lines are interchanged (so that the external lines enter in t)' and exit from 0, see Figure 
6) the momentum through the line £i changes from i/ e o + v to —vgs + v. 




Figure 6. The self-energy cluster D' obtained from D by interchanging the external 
lines of D: the exiting line is attached to a node € D' 2 and the entering line is 
attached to a node ti£DJ. 

At x = 0, if we change D into D' , the value V£) 1 (0;?7)9 x g , [™ < il(a;- ) ;r/)V£) 2 (0;ry) is changed into 
Vd'(0;?7) d x g^- ne ^ (— x®; 7])Vd' (0; if). In analogy with the results of Appendix A5, employing the 
symmetry relations (A5.1), we find V d > 2 (0;t]) = E[Vd 2 (0; r])] T E and Vo^rf) = E[V Dl (0;n)] T E. 
Using also the symmetry Ed x g^ ne ^(—x^;r])E = —d x [g^ ne ^] T (x^; rj), we find that the element 
(Ai, otj) of the contribution corresponding to the graph in Figure 6 can be written as 

(v D;2 (0; V )d x g^(~x^; V )V D[ (0; V )) = -(E[V D2 (0;v)} T d x [g [n ^} T (^v)Vl 1 (^v)E) a = 

= -(EV Dl (0;v)d x g ln ^(x l ;r 1 )V D2 (0;r 1 )E) = -(v Dl (0;v)d x g [n ^(^;v)VD 2 (0;v)) A ■ 

(A6.5) 

Repeating the discussion leading to (A6.4) we see that, if we sum (A6.5) over the graphs obtained 
by attaching the entering line to all possible nodes t)' E D[ Ci D and attaching the exiting line to 
all possible nodes t) € D' 2 n D , we get 

-( E M 3 )^M,dd x 9 ln ^\^;v))^M^% 2 ( E 046.6) 

Now, using the fact that X)oen ^0=0 and that D[,D' 2 and D\,D 2 are topologically equivalent, 
we also find that 

E { v *')i = - E E (^)* = - E ( A6J ) 

o'eDJnflo 0G-D2n_D o 0G-D 2 nr> o'eDinDo 
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so that (A6.4) and (A6.6) sum up to 0. 
This completes the proof of the fact that 

V d(£: 0)A iaj = 0, 8 X Vd(e; x) Aiaj | x=Q = 0. (A6.8) 

def 

Remark. It can be checked that the symmetry (A5.1) "only" implies that, if we define Cij(x) = 
^Der Vd(£', x )A iaj , it holds Cij(x) = Cji(-x) = C*j(x): therefore the cancellation to second order 
(in x), expressed by (A6.7), is not a parity cancellation (unless i = j where by the self-adjointness 
of the matrix C it follows that Cu(x) is real and, hence, even in x). If the coefficients f v (f3) are 
assumed real then the self-adjointness property (A5.1) implies that Cij(x) is even in x and the 
second order cancellation at x = is an obvious parity cancellation and the proof above is greatly 
simplified, as remarked in the simple cases considered in [Gl] and [G2]. 

Having proved that (M [n] (0; rj)) , = if either 7' = A; or 7 = aj and that (d x M [n] (0; 77)) , = 
if 7' = Ai and 7 = ctj, it has still to be shown that the elements (M^ (x; f?)) , (with the 
suitable choices of 7' and 7) satisfy appropriate bounds once the factors x determining the order 
of zero at x — are extracted. From convergence one expects that the bounds on the rests of the 
Taylor series in x around x = should still be proportional to e 2 . 

This is in fact true, under the condition that (A6.1) does not hold: if (A6.1) does not hold, then 
changing the nodes where the external lines are attached to does not change the scale of the internal 
lines. And this implies that the bounds on the derivatives with respect to x are qualitatively the 
same as the bounds on M\ n ^(x\ if) and boundedness of the entries of ikf'~"' follows, see (5.12). 

Appendix A7. Bounds on the propagator 

In this section we want to get a lower bound for the norm of (ix — M M (x; rj))E, by an approximate 
computation of its eigenvalue which is lowest in absolute value, for x 2 < pry 2 * 0-1 , with p the same 
constant appearing in Lemma 3. Moreover we want to prove dimensional bounds for the derivatives 
in x and e of the approximate lowest eigenvalue. 

Write {ix M^(x; V ))E as (ix M^(x; V ))E = [ ix ~+Q[<^ v) '^piM^'^ 

where pl< n \ Ql< n \ R[< n 1 are the matrices defined in (5.12). 

The N eigenvalues of P^- n \x;rj) are of order 1 and positive; call them fj,^ < ... < pff. 
Since P^- n \x;r]) differs from P^(x;n) by 0(e 2 ), the eigenvalues can be written as fj,^ = 
M, (0) (l + O(e 2 )). 

To estimate the smallest eigenvalue of (ix — M\- n \x;ri))E we can follow a strategy adapted 
from the proof of Lemma 2. Let v a (x;rj) — ( U ) be an eigenvector of (ix — M\- n \x;rj))E with 



(47.1) 

u 



eigenvalue a (here u and v are two column vectors of dimension N): 



-R^(x;v) -ix + Q^ n \x;r)) 
ix + Q^ n ^(x-r 1 ) Pte"\{ X ;rj) 

_ ( -R^u + (-ix + Q^)v 
(ix + Q^ n ^)u + P^ n ' 



v 



a 

VI it) 
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To estimate the eigenvalue a with smallest absolute value, we restrict attention to the case 

p[<n] _ 

a > p[<"]/2, as in the opposite case |a| > m^V 2 - Thcn {P [ - n] - a) is invertible, (P^"l - a)- 1 = 
[p[<™]]-! -). O(a) and we can rewrite (A7.1) as 



- P^"l - (-is + QK"])(P^"] - a)" 1 (is + Q^)u = au. 
From the second of (A7.2) we see that a must solve an equation of the form 



(A7.2) 



dct 



- a(l + O(x) + 0(e))] = 0, (A7.3) 



where we defined d = f -R^ - (-ix + Q^^P^-^x + Q^t) and we used that Q^"l = 

0(e). This means that a is an eigenvalue of a Hermitian matrix of the form JV^-"! (l + 0(x) + 0(e)) , 
so that, calling A'"' the eigenvalues of N^- n \ we have (see [Ka]) 

a= X [ ; i] (l + 0(x)+0(e)), (A7A) 

for some i = 1, . . . , N. The conclusion of the previous discussion is that the problem of computing 
the smallest eigenvalue of (ix — (x; rj))E is essentially equivalent to the problem of computing 
the smallest eigenvalue of iV'- n '. 
An explicit computation of N^- n ^ shows that it can be written in the form 

N [<n] = _ X 2 H -1 + | _ R [<n] + X 2 H -1 _ ^ + Q [<n] ) [p[<„] ] - 1 ^ + g[<«]t)} = 

where W ( ) , %(,,) = -cfco/jf""^! + Ofa)) and D( V ) = D + 0( V ), with 

The off-diagonal O(ex) elements can be eliminated through a small rotation of N^- n h In fact, 
defining the Hermitian matrix Z as: 



z= ( _ ^(i,)]- 1 ^) 

{[M^ri)]-^) o 

a computation shows that 

elxv -^Z N [<n] e - lx11 -^Z = _a.2 w -l + 5JV [<n] = 



(A7.7) 
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where we used the property xr/~ ko+1 < (pry) 1 / 2 , if x 2 < perj k °^ 1 . Multiplying the r.h.s. of (A7.8) 
times Ti.^ 2 both from the left and from the right (remark that H^j 2 is a well-defined positive 0(1) 
operator) we see that the norm of jVl<«] 

can be bounded above and below by an 0(1) constant 

times the norm of 

-x 2 +nl /2 SN^nl /2 . (A7.9) 

Finally using the explicit structure of 5N^- n \ see (A7.8), N — 1 of the eigenvalues of (A7.9) are 
equal to -x 2 (l + Oil] 1 / 2 )), while the last one is -x 2 (l + 0(r/ 1/2 )) - en ko - 1 Mf 3 f 3 (s)dlH Q . 

From this computation, the bound (6.1) of Lemma 5 follows. To address the differentiability 
of \W(x;ri) with respect to e and x, we note that H^ 2 6N^- n ^'H^ 2 is C°° in the sense of Whit- 
ney (because such is M\- n \ see the corresponding discussion in Appendix A3 of [GG2]) and its 
derivatives admit the following dimensional bounds: 

\\d e (nl /2 SN^nl /2 )\\ <c v ko - 1 , \\d x (nl /2 SN^nl /2 )\\ <c v ko . (alio) 

Given this and the fact that the last eigenvalue of H^ 2 SN^-'^H^ 2 is isolated (it is non degenerate 
and its distance from the others is O(eri ko ^ 1 )), we can represent X^(x; if) in the form: 

A™ (x; rj) = -rf^^Tr f I ^ ^] , (A7.ll) 

where z is a complex variable independent of e, x and 7 is a circle in the complex plane around 
Mpp^dgHo not sorrounding the origin, with e-independent radius. The derivatives of \^ n \x;rj) 
with respect to e and x can be computed differentiating the r.h.s. of (A7.ll) and, using the 
dimensional bounds (A7.10), we find the same dimensional bounds for the derivatives of A^ (a;; rj): 

\d E \ ln] (x;n)\ < Cr) k °-\ \d x \ [n] (x; rj)\ < Cn ko . (A7.12) 



The proof of Lemma 5 is concluded. 



Appendix A8. Excluded values of the perturbation parameter 

In this section we want to check that, by imposing the first condition in (6.8), the measure of 
excluded values of e can be bounded by the second of (6.8). To estimate the measure of the 
excluded e's we proceed as in Appendix A2 of [GG2]. We present a proof valid for any n > n^. 
The dimensional bounds (6.2) and the property |A^(7j) - At" - 1 ]^) < Ce~ Kl2 " /T1 (see Lemma 5) 
implies the following dimensional bound for the derivative of \J X^(r]): 



\d e Jmv)\<c'r y2 - 048.i) 

Then the first condition in (6.8) excludes, for each is, a subinterval of 1(e) whose measure is 
bounded by 

(~i 9-m/2„l/2 
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The Diophantinc condition on uj implies that if the first condition in (6.8) is invalid then \v\ cannot 
be too small: 

+ > \x\ > CoW\- T » . (A8.3) 

Therefore > CV°~ 1/2 , hence in this case we only have to consider the values of u such 

that \v\ T ° > (Co/2C")ri^~ ko . Summing (A8.2) over the z/s satisfying this constraint, and using 
the definition n = r (l + (5i) + N we find that the total measure of excluded e's can be bounded 

which is in fact the second inequality in (6.8). 
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